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Abstract 

We calculate a renormalized Hamiltonian for pure-glue QCD and diagonalize it. The renor- 
malization procedure is designed to produce a Hamiltonian that will yield physical states that 
rapidly converge in an expansion in free-particle Fock-space sectors. To make this possible, 
we use light-front field theory to isolate vacuum effects, and we place a smooth cutoff on the 
Hamiltonian to force its free-state matrix elements to quickly decrease as the difference of 
the free masses of the states increases. The cutoff violates a number of physical principles of 
light-front pure-glue QCD, including Lorentz covariance and gauge covariance. This means 
that the operators in the Hamiltonian are not required to respect these physical principles. 
However, by requiring the Hamiltonian to produce cutoff-independent physical quantities and 
by requiring it to respect the unviolated physical principles of pure-glue QCD, we are able to 
derive recursion relations that define the Hamiltonian to all orders in perturbation theory in 
terms of the running coupling. We approximate all physical states as two-gluon states (thus 
they are relatively simple single-glueball states), and use our recursion relations to calculate 
to second order the part of the Hamiltonian that is required to compute the spectrum. We 
diagonalize the Hamiltonian using basis-function expansions for the gluons' color, spin, and 
momentum degrees of freedom. We examine the sensitivity of our results to the cutoff and 
use them to analyze the nonperturbative scale dependence of the coupling. We investigate 
the effect of the dynamical rotational symmetry of light-front field theory on the rotational 
degeneracies of the spectrum and compare the spectrum to recent lattice results. Finally, we 
examine our wave functions and analyze the various sources of error in our calculation. 
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1 Introduction 



A solution to a quantum field theory that is close to our physical intuition is possible if we can develop 
a formalism in which the physical states of the theory rapidly converge in an expansion in free-particle 
Fock-space sectors (free sectors). Such an expansion is unlikely to be possible in an equal-time approach 
to many of the more interesting theories, such as QCD. This is because in equal-time approaches to these 
theories, the physical states must be built on top of a complicated vacuum unless the volume of space 
is severely limited. For this reason, we work in light-front field theory (LFFT). In LFFT, it is possible 
to force the vacuum to be empty by removing from the theory all particles that have zero longitudinal 
momentum|^. Any physical effects of these particles must be incorporated into the operators of the theory 
in order to obtain correct physical quantities. 

In LFFT, the Hamiltonian is trivially related to the invariant-mass operator^ (IMO), and it is more 
natural to work with the IMO because it is manifestly boost-invariant. If the IMO satisfies three conditions 
in the basis of free-particle Fock-space states (free states), then its eigenstates will rapidly converge in 
an expansion in free sectors. First, the diagonal matrix elements of the IMO must be dominated by the 
free part of the IMO. Second, the off-diagonal matrix elements of the IMO must quickly decrease as the 
difference of the free masses of the states increases. If the IMO satisfies these first two conditions, then 
each of its eigenstates will be dominated by free-state components with free masses that are close to 
the mass of the eigenstate. The third condition on the IMO is that the free mass of a free state must 
quickly increase as the number of particles in the state increases. If the IMO satisfies all three conditions, 
then the number of particles in a free-state component that dominates an eigenstate will be limited from 
above. This means that the IMO's eigenstates will rapidly converge in an expansion in free sectors^. 

To satisfy the first condition on the IMO, we assume that we can derive the IMO in perturbation 
theory. If this is true, then the couplings are small and the diagonal matrix elements of the IMO are 
dominated by the free part of the IMO. To satisfy the second condition, we place a smooth cutoff on the 
IMO to force its matrix elements to quickly decrease as the difference of the free masses of the states 
increases. Once we have removed the particles with zero longitudinal momentum from the theory, it is 

^This is because there are no negative longitudinal momenta and momentum conservation requires the three- momenta 
of the constituents of the vacuum to sum to zero. 

^The invariant-mass operator is given by the square of the momentum operator: V^V^ = M^. See Appendix A for 
more details. 

^There are three subtleties here. The first subtlety is that the first and third conditions on the IMO will not be satisfied 
for those free states in which many of the particles have negligible center-of-mass transverse momentum and little or 
no mass. However, the contributions of these free states to the physical states in which we are interested are typically 
suppressed. For example, in QCD these free states have very large widths in transverse position space and are thus highly 
suppressed by confinement. In QED, the particles with negligible center-of-mass transverse momentum and no mass are 
long-transverse-wavelength photons. These photons decouple from the physical states in which we are typically interested, 
e.g. charge-singlet states like hydrogen and positronium. Thus the contributions to these physical states from the free states 
containing these photons are suppressed. 

The second subtlety is that exactly how quickly the IMO's off-diagonal matrix elements must decrease and the free mass 
of a free state must increase are not known. We assume that the rates that we are able to achieve are sufficient. This can 
be verified by diagonalizing the IMO and examining the rate of convergence of the free-sector expansion of its eigenstates. 

The third subtlety is that the coefficients of the expansion for highly excited eigenstates may grow for a number of free 
sectors and then peak before diminishing and becoming rapidly convergent. 
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reasonable to expect that the third and final condition on the IMO will be satisfied automatically. This is 
because the free-particle dispersion relation of LFFT should force the free mass of a free state to quickly 
increase as the number of particles in the state increases (see Appendix A of Ref . ) . 

By suppressing the matrix elements of the IMO that have large changes in free mass, the cutoff 
regulates the ultraviolet divergences of the theory. Unfortunately, it also violates a number of physical 
principles of LFFT, including Lorentz covariance and gauge covariance. This means that the operators 
in the IMO are not required to respect these physical principles, and renormalization is not simply a 
matter of adjusting a few canonical parameters. The simplest way to systematically determine the IMO 
in this case is in perturbation theory. In order for a pcrturbative computation of the IMO to be strictly 
valid, the theory must be asymptotically freeQ. If this is the case, then by requiring the IMO to produce 
cutoff-independent physical quantities and by requiring it to respect the unviolated physical principles of 
the theory^, we can derive recursion relations that define the IMO to all orders in perturbation theory in 
terms of the fundamental parameters of the field theory. If our cutoff is large enough, then the couplings 
will be small and the pcrturbative approximation to the IMO may work well. 

The physical principles that we use to determine the IMO form a subset of the full set of physical 
principles of light-front field theory. This raises the question of how the remaining principles, which are 
violated by our cutoff, are restored in physical quantities. Since the IMO is uniquely determined by the 
principles that we use, the remaining principles must be automatically respected by physical quantities 
derived from our IMO, at least perturbatively. If they are not, then they contradict the principles that 
we use and no consistent theory can be built upon the complete set of principles. The reason that this 
process is possible is that there are redundancies among the various physical principles. 

It is possible to compute operators other than the IMO in our approach. Although wc compute op- 
erators perturbatively, we can use these operators to compute nonperturbative quantities. For example, 
the spectrum can be computed by diagonalizing the IMO (see Section 6). However, there are drawbacks 
to computing operators perturbatively. It is possible that there are intrinsically nonperturbative effects 
in the theory that require nonperturbative renormalization. Any such effects are neglected in this ap- 
proach. Another problem is that perturbative renormalization makes nonperturbative physical quantities 
somewhat cutoff-dependent. 

In general, field theories have an infinite number of degrees of freedom. However, since our IMO 
will cause the physical states of the theory to rapidly converge in an expansion in free sectors, we can 
truncate this expansion. This means that approximate computations of physical quantities will require 
only a finite number of finite-body matrix elements of operators. In addition, since we assume that we can 
compute these matrix elements perturbatively, we do not implement particle creation and annihilation 
nonperturbatively in these calculations. This allows us to always work with a finite number of degrees of 
freedom. 

When we perturbatively calculate matrix elements of operators, we do not truncate the space of 

*Our method may work even if the theory is not asymptotically free. For example, it works in QED because the scale 
at which the electron charge is large is astronomical. 

^Some of the physical principles, such as cluster decomposition, are violated in a very specific manner and can still be 
used to restrict the form of the IMO. However, the restriction in this situation is always weaker than it would have been 
had there been no violation of the principle. 
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intermediate states that can appear; so our approach does not use a Tamm-Dancoff truncation^ ^, 
We also do not completely eliminate any interactions, such as those that change particle number. These 
strengths of our approach allow us to better describe physical theories. However, the truncation of the free- 
sector expansion of physical states has drawbacks that are similar to those of perturbative renormalization. 
It neglects any physical effects that require an infinite number of particles and contributes to the cutoff 
dependence of nonperturbative physical quantities. 

The accuracy of our results and the strength of the cutoff dependence of our nonperturbative physical 
quantities are determined by the order in perturbation theory to which we calculate the operators of the 
theory and the number of free sectors that we keep in the expansion of physical states. If we use a cutoff 
that is too small, then the couplings of the theory will be large, and it will be necessary for us to keep 
many terms in the expansion of the operators. If we use a cutoff that is too large, then the free-sector 
expansion of the physical states will converge slowly, and it will be necessary for us to keep many sectors 
in the states. We assume that if the order of perturbation theory and the number of free sectors are 
manageable, then there is a range of cutoff values for which the approximations work well and physical 
quantities are relatively accurate and cutoff- independent. 

As we mentioned, we remove from the theory all particles with zero longitudinal momentum. We 
should replace their physical effects with interactions. However, due to the limitations of our method, we 
can reproduce only those effects of these particles that can be derived with perturbative renormalization 
and require only a small number of particles. 

There are a number of approaches that are similar to ours j|, ^, 0, |[ |9[ |ll|, |l^ , and some of these 
methods have been used to calculate the physical states of QCD These calculations are based 

on nonrelativistic approximations and use sharp step-function cutoffs. Nonrelativistic approximations 
drastically simplify the diagonalization of the IMO, but are insufficient for states containing light quarks 
or giuons. Sharp cutoffs prevent the complete cancellation of the infrared divergences that appear in 
light-front gauge theories^]. 

Our approach is completely relativistic and uses smooth cutoffs to ensure the complete cancellation 
of the light-front infrared divergences. It is largely based on the renormalization methods of Perry |p^ . 
Perry and Wilson ||l^ , Wilson , and Glazek and Wilson ||^ , as well as the Hamiltonian-diagonalization 
methods of Wegner |17 . In Ref. Q], we developed and tested our method in massless (fP' theory in six 
dimensions. In Ref. jl^, Kylin, Allen, and Perry extended this method to include particle masses. In this 
paper, we extend our method to pure-glue QCD. This theory is simpler than full QCD due to the reduced 
number of vertices and the absence of quark masses. For this application, we derive the recursion relations 
that determine the IMO to all orders in perturbation theory in the running coupling. We approximate all 
physical states as two-gluon states and use the recursion relations to compute to second order the part of 
the IMO that is required to compute the spectrum. We diagonalize the IMO in a basis-state expansion 
and analyze the results. 

®The only degrees of freedom that we remove from the full theory are the particles with zero longitudinal momentum. 
We should be able to replace the physical effects of these particles with interactions without compromising the validity of 
the theory. This is because these particles are vacuum effects or have infinite kinetic energies (or both), and thus are not 
observable in the laboratory as particles. 

''These appear due to the exchange of massless gauge particles with arbitrarily small longitudinal momentum. 
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In the remainder of the Introduction, we outline the rest of the paper. Appendix A contains our 
conventions for light-front pure-glue QCD; so the reader may wish to examine it first. In Section 2 we 
present our method for computing the free-state matrix elements of the IMO. After defining the action 
of our cutoff, we place a number of restrictions on the IMO. First wc force the IMO at a given cutoff 
to be unitarily equivalent to itself at a higher cutoff. This implies that the IMO is unitarily equivalent 
to itself at an infinite cutoff, and will therefore yield cutoff-independent physical quantities. From the 
statement of unitary equivalence, we develop a perturbative series that relates the interactions at two 
different cutoffs. We then proceed to use physical principles to restrict the form of the IMO. We require 
it to conserve momentum and to be invariant under boosts and rotations about the 3-axis. Although 
our cutoff violates exact transverse locality, we arc able to require the IMO to respect an approximate 
transverse locality. In practice this means that the IMO's matrix elements arc analytic functions of 
transverse momenta. The cutoff also violates cluster decomposition, but we show that the implications 
of this violation are simple enough that we can still use this principle to restrict the form of the IMO. 

In order to represent pure-glue QCD, the IMO must become the free IMO in the noninteracting limit, 
and it must be a function only of the cutoff and the coupling. The final physical restriction on the IMO is 
that it must reproduce the perturbative scattering amplitudes of pure-glue QCD. This restriction specifies 
the form of the first-order interaction and part of the second-order interaction. 

We use the restrictions from physical principles and the perturbative series that relates the interactions 
at two different cutoffs to derive the recursion relations that determine the IMO. A crucial step in this 
process is the removal of the coupling from the perturbative series. This allows us to separate the cutoff 
dependences of the operators in the interaction from the cutoff dependences of their couplings. Most of 
the details of this derivation are relegated to Appendix B. 

In Section 3 we define a basis for the expansion of physical states. We construct the physical states 
to be simultaneous eigenstatcs of the IMO, the three-momentum operator, and the part of the generator 
of rotations about the 3-axis that generates rotations of the internal degrees of freedom of states (as 
opposed to their centers of mass). We approximate all physical states as two-gluon states, which means 
that they are relatively simple single-glueball states. We analytically calculate the color parts of the 
states by requiring the states to be color singlets, and then we expand the states' spin and momentum 
degrees of freedom in a complete set of orthonormal basis functions. We conclude Section 3 by deriving 
the eigenvalue equation for the IMO in our basis. 

In Section 4 we use the recursion relations for the IMO to compute its two-gluon to two-gluon matrix 
element to second order in the running coupling. This is the only free-state matrix element that we need 
to solve the eigenvalue equation. 

In Section 5 we compute the matrix elements of the IMO in the basis that we have defined, in terms 
of integrals that must be evaluated numerically. To avoid roundoff error, we evaluate the integrals for the 
kinetic energy and the two-point interaction (the self-energy) by writing them as sums of gamma functions. 
The remaining integrals are five-dimensional, and we compute them with Monte Carlo methods. However, 
before we can do this, it is necessary for us to make manifest the cancellation of the infrared divergences 
from exchanged gluons with infinitesimal longitudinal momentum. This is nontrivial, and even after 
the divergences are gone, we must manipulate the integrals quite a bit to get them into a form that is 
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amenable to Monte Carlo integration. We present some of these details and other technical issues in 
Appendix C. 

In Section 6 we derive our results and analyze the sources of error in the calculation. We begin by 
discussing how we assign quantum numbers to our numerical results for glueball states and proceed with 
a discussion of the procedure that we use to compute these results. We derive the nonperturbative cutoff 
dependence of the coupling and discuss the cutoff dependence of our glueball masses. We use this analysis 
to choose the value of the cutoff that minimizes our errors. We then present the spectrum that we find 
with this optimal cutoff and compare it to recent lattice results. The last results that we present are the 
probability densities for our five lightest glueballs. We conclude Section 6 by discussing the sources of 
error in our calculation and estimating the sizes of these errors. 

Finally, in Section 7 we conclude with a summary and a discussion of the direction of future work. 

2 The Method for Computing Free-State Matrix Elements of 
the Invariant-Mass Operator 

This section summarizes the extension of the results of Sections 2-4 and Appendix D of Ref. to the 
case of pure-glue QCD. 

2.1 The Cutoff 

Our goal is to derive recursion relations that uniquely determine the IMO to all orders in perturbation 
theory in the running coupling. The IMO is a function of the cutoff, and can be split into the canonical 
free IMO and an interaction (see Appendix A for our light-front pure-glue QCD conventions): 

M\A)^ML, + MUA). (1) 

The cutoff is implemented on the matrix elements of the IMO: 

{F\M'iA)\I) = {F\MLJI) + {F\MUA)\I) 

= Mj{F\I)+e-^{F\ViA)\I), (2) 

where \F) and |/) are eigenstates of the free IMO with eigenvalues Mp and Mf, and Apj is the difference 
of these eigenvalues: 

Apj = Ml-Mf. (3) 

V{A) is the interaction with the Gaussian cutoff factor removed, and we refer to it as the "reduced 
interaction." To determine the IMO, we must determine the reduced interaction. 

We will see that (^F\ V^(A) |/) does not grow exponentially as A'pj gets larger so the exponential in Eq. 
forces the off-diagonal matrix elements of the IMO to rapidly diminish as Aj^j grows. This satisfies 
the second of our three conditions on the IMO and regulates it. 
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2.2 The Restriction to Produce Cutoff-Independent Physical Quantities 

Our cutoff violates a number of physical principles of LFFT, including Lorentz covariance and gauge 
covariancc^. This means that the operators in the IMO are not required to respect these physical 
principles. In addition, since there is no locality in the longitudinal direction in Hamiltonian LFFT|^ 
these operators can contain arbitrary functions of longitudinal momenta. To uniquely determine the 
IMO in this case, we have to place some restrictions on it. The first restriction is that it must produce 
cutoff-independent physical quantities. To enforce this, we require A4^(A) to satisfy 

M^{K) = U{K,K')M^{h')U\K,M), (4) 

where U is a. unitary transformation that changes the IMO's cutoff and 

A < A' < 2A. (5) 

We have placed an upper limit on A' because Eq. (Q) is perturbatively valid only if A' is not too much 
larger than A [|l9j. Note that we are considering A^'^(A) to be a function of its argument; i.e. A^^(A') 
has the same functional dependence on A' that A^^(A) has on A. In Rcf. we proved that Eq. (||) 
forces A^^(A) to produce cutoff- independent physical quantities. 

The unitary transformation that we use is designed to alter the cutoff implemented in Eq. (|^), and 
is a simplified version of a transformation introduced by Wegner [ pTf , modified for implementation with 
the IMO. It is uniquely defined by a linear first-order differential equation: 

dC/(A,A') 



d(A-4) 



T(A){/(A,A'), (6) 



with one boundary condition: 



[/(A,A) = 1. (7) 
f7(A, A') is unitary as long as r(A) is anti-Hermitian and linear We define 

T(A)=[A^L„A<2(A)], (8) 

which is anti-Hcrmitian and linear. 

To solve for A^^(A) perturbatively, we need to turn Eq. (||) into a perturbative restriction on the 
reduced interaction, V^(A). We outlined how to do this in Ref. and here we simply state the results 
that we need. The perturbative version of Eq. (Q) in terms of the reduced interaction is 

V{K)^V{K') = 5V, (9) 
where 5V is the change to the reduced interaction and is a function of both A and A': 

®Our regulator breaks these symmetries because the mass of a free state is neither gauge-invariant nor rotationally 
invariant (except for rotations about the 3-axis). 

"That there is no longitudinal locality in Hamiltonian LFFT is evident from the fact that the longitudinal momentum 
of a free particle appears in the denominator of its dispersion relation, (joj + m^)/p+. 
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{F\SV\I) ^ i5](^^|nA')|i^)(i^inA')|/)rf''''^(F,if,/) 

K 

+ i ^ (F| ^(A') \K) {K\ V{A') \L) {L\ F(A') |/) Tf ^^')(^^, K, L, I) 

K,L 

+ o([ViA')t). (10) 



In this equation, the sums are over complete sets of free states and the cutoff functions are defined by 

1 1 



and 
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1 1 



Aif/ Ap^ 



g2A'-''A^i,Ai, ^g2A-»Ay^Aj^, _ g2A'-''Aj,^ A^, ~j 



^FK + A 



FA- ^Jg |^g2A'-''(A^^.A^,+A^^A^,) _ g2A-^(Ap^A^,+A,,.^A^,) 



\^FK ^KlJ \^LI A^J V J 



AfL + A 
A^^A,^, + A^,A 



FL ^IL |^g2A'-''(Ap^A^^+A^^A^,) _ g2A-^(A^^ A^-^+A^^A^,) 



(12) 



The above definitions for the cutoff functions assume that none of the A 's that appear in the denominators 
is zero. In the event one of them is zero, the appropriate cutoff function is defined by 

r/A^^')(A=0) = Hm^r/^''^'^(A). (13) 
2.3 Restrictions from Physical Principles 

Eq. (^) is the first restriction on the IMO. To imiquely determine the IMO, we need to place additional 
restrictions on it, and we do this using the physical principles of the theory that are not violated by the 
cutoff. (See Footnote ||.) 

2.3.1 Symmetry Principles 

Any LFFT should exhibit manifest momentum conservation, boost covariance, and covariance under 
rotations about the 3-axis. Our cutoff does not violate any of these principles; so we restrict the IMO to 
conserve momentum and to be invariant under boosts and rotations about the 3-axis. 



2.3.2 Transverse Locality 

Ideally, the IMO should be local in the transverse directions, and thus each of its matrix elements 
should be expressible as a finite series of powers of transverse momenta with expansion coefficients 
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that arc functions of longitudinal momenta. In our case, the cutoff suppresses interactions that have 
large transverse-momentum transfers and replaces them with interactions that have smaller transverse- 
momentum transfers. This is equivalent to suppressing interactions that occur over small transverse 
separations and replacing them with interactions that occur over larger transverse separations; so we 
do not expect our interactions to be perfectly transverse-local. Nonetheless, we expect that interactions 
in A^^(A) should appear local relative to transverse separations larger than or, equivalently, to 
transverse momenta less than A. This means that for transverse momenta less than A we should be 
able to approximate each matrix element of A^^(A) as a finite power series in p^/A. We enforce this by 
assuming that transverse locality is violated in the weakest manner possible, i.e. that any matrix element 
of the IMO can be expressed as an infinite series of powers of transverse momenta with an infinite radius 
of convergence. In other words, we assume that the matrix elements of the IMO are analytic functions 
of transverse momenta. 

2.3.3 Cluster Decomposition 

Since the matrix elements of the IMO conserve momentum, they can be written as a sum of terms, with 
each term containing a unique product of momentum-conserving delta functions |^o| . We require the IMO 
to satisfy approximate cluster decomposition Q; i.e. when any of its matrix elements is written as an 
expansion in the possible products of momentum-conserving delta functions, the coefficient of any term 
in the expansion is restricted as follows. (If there is more than one possible set of spectators for a given 
product of momentum-conserving delta functions, then the coefficient has to be broken into a distinct 
part for each possible set, and these restrictions hold for each part separately.) It can depend on the 
cutoff, the quantum numbers of the interacting particles, and the total longitudinal momentum. It must 
be proportional to a quantum-number-conserving Kronecker delta for each discrete quantum number for 
each spectator. It can have no other dependence on the quantum numbers of spectators, and it cannot 
contain delta functions that fix momenta. 

The reason that the IMO does not respect exact cluster decomposition and that the coefficients in 
the delta-function expansion can depend on the momenta of spectators (through a dependence on the 
total longitudinal momentum) is that our cutoff on free-mass differences violates cluster decomposition. 
To see this, note that the change in free mass for some process is given by P^{J2jP'j^ ~ J2iP7)j where 
the Pi's are the momenta of the particles in the initial state, the p^s are the momenta of the particles 
in the final state, and P+ is the total longitudinal momentum of each state. The minus momenta of any 
spectators cancel in this difference, but their longitudinal momenta still contribute to the overall factor 
of P+. 

2.3.4 Representation of the Theory of Interest 

The preceding restrictions on A^^(A) are valid for any LFFT in more than two dimensions. In order to 
represent a particular theory, we must place additional restrictions on M'^(A). 

We assume that we can compute A^^(A) perturbatively, which means that we can expand T^(A) in 
powers of the coupling at the scale A. Our cutoff has no effect in the noninteracting limit; so our IMO 
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must reproduce free pure-glue QCD in this limit. According to Eq. (H), this means that ^(A) vanishes 
in the noninteracting limit. 

In pure-glue QCD, the only fundamental parameter is the coupling; so we require the IMO to depend 
only on it and the scale. (For an example of the application of our method to a theory with more than 
one parameter, see Ref. [fsl.) In this case, the expansion of V{A) takes the form 



nA) = £5:^W(A), 



(14) 



r=l 



where is the coupling at the scale A. We refer to ^(''^(A) as the 0{g^) reduced interaction, although 
for convenience the coupling is factored out. 

is the correct fundamental parameter for pure-glue QCD if and only if its definition is consistent 
with the canonical definition of the coupling. The canonical definition is 

-1 
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167r3p+5(3)(p^_P2-p3) -^^9293^91^ {9293\ \gi) , 



(15) 



where v is the canonical interaction and <|;j|u|z^ denotes a modified matrix element of v. (See Appendix A 
for our light- front pure-glue QCD conventions.) The denominator removes all dependence on momentum, 
spin, and color in the canonical matrix element for giuon emission, and thus isolates the coupling. We 
define the coupling by 

-1 ~ -1 . - 
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167r3p+5(3)(p^ „P2 - Pa) <92.93|w|gi:t> 

i6n^ptS^^\pi-P2-P3) ^g2g3\v\gi:^ 



(32.931 A^'(A) |5i 



Sti— 1; Cn—n 



(16) 



P2±=P3±; pj=pj ; e=o 

where n — l,2,S, and e is a longitudinal- momentum cutoff that we define below. Momentum conservation 
and boost invariance imply that the matrix element that defines the coupling can depend only on the 
transverse momentum of particle 2 in the center-of-mass frame and the ratio P2 /Pi ■ The restrictions 
P2i_ = P3i_ and P2 = p^ fix these quantities to be and 1/2, respectively. 

Our definition of the coupling is consistent with the canonical definition because the conditions on 
the matrix elements in Eq. ( |l6|) have no effect on the right-hand side (RHS) of Eq. ( [l5| ) and do not force 
(5253!-^^ to vanish. According to Eq. ([l^), the IMO is coupling coherent [|l^, |l^, |l| because the 
couplings of its noncanonical operators are functions only of the fundamental parameters of the theory 
and they vanish in the noninteracting limit. 

We have assumed that the IMO obeys approximate transverse locality, which means that we can 
expand any matrix element {F\ ^'^''''(A) |/) in powers of transverse momenta. Each term in this expansion 
is either cutoff-dependent or cutoff-independent. We define V^^{h.) to be the cutoff-dependent part of 
y('")(A), i.e. the part that produces the cutoff-dependent terms in transverse-momentum expansions of 
matrix elements of ^(''^(A). We define V^j'* to be the cutoff-independent part of y^''^(A), i.e. the part 
that produces the cutoff-independent terms in transverse-momentum expansions of matrix elements of 
1/('~)(A). Then 



y«(A) = v^W(A) + yW. 



(17) 
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This separation is necessary because the procedures for computing V'qp(A) and Vq^ differ. 

We beheve that if A^^(A) is to reproduce the perturbative scattering amphtudes of pure-glue QCD, 
then it is necessary and sufficient for the reduced interaction to satisfy the following condition^^. We 
prohibit ^('"^(A) from having a three-point interaction unless r is odd, and we prohibit y('')(A) fr om 
having a four-point interaction unless r is even. We require 

^(1) ^ j D^D2D:ialalai5^^\p^-p2-pz)^g2gM9i~k 

+ P2 -pa) "i;5'3|w|5'i52^ (18) 



and 



^CI 



where 



and 



and 



■ ■ J i=l 

'P^^r / ^1^2£'3£'4 444^1 <5'^Hp1 -P2 -P3 -P4) 51 ^ll3,4 ^929:i9i\v\9l> ^ 
■ i=l 

■ / DiD2D3D4alaia2a3S^^^(pi+p2+P3-P4)^0il^3.4 ^941^1919293^1, (19) 



Wit- 



,(1) 

1,2:3,4 



1, 



^tlsA = 0{\pt-pt\-er+), 
5(3) - ati^+_^+ 



il3,4 = 0{\pt~pt\-eV+), 

1, (20) 



,(4) 
1,2:3,4 



,(1) 
1;2,3,4 

9(2) 
1;2,3,4 

9(3) 
1;2,3,4 

9(4) 
1;2,3,4 



9(1) 

1,2,3;4 

9(2) 
1,2,3;4 

9(3) 
1,2,3;4 

9(4) 
1,2,3;4 



1, 






0{\pt 


-pt\ 




o{\pt 


-Pt\ 


- er+), 


Oi\pt 


-pt\ 





(21) 



1, 






0{\pt 


-pt\ 




o{\pt 


~pt\ 


- eV+), 


0{\pt 


-pt\ 


-eV+). 



(22) 



The presence of e in these step-function cutoffs ensures that we will avoid divergences from exchanged 
gluons (either instantaneous or real) with infinitesimal longitudinal momentum. In Section ^ we show 
how we can take e to zero before we diagonalize M'^{A). 

^"Any longitudinal regulator that is consistent with the physical principles that we use to restrict the IMO is sufficient. 
It is not necessary to use our e cutoff in order to reproduce the perturbative scattering amplitudes of pure-glue QCD. 
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We have not yet proved that the above conditions on the reduced interaction are necessary and 
sufficient to reproduce the perturbative scattering aniphtudes of pure-glue QCD. We can, however, show 
in (jr' theory 0] that analogous conditions are necessary and sufficient to reproduce all second-order 
scattering amplitudes. 

2.4 The Recursion Relations for the Invariant-Mass Operator 

The restrictions that we have placed on the IMO are sufficient to uniquely determine it order-by-order in 
perturbation theory. In this subsection, we present the recursion relations that define A^^(A) in terms of 
and V^^, which we have defined above. To begin, we consider the restriction that forces the IMO 
to produce cutoff-independent physical quantities: 

y(A)-y(A') = 5V. (23) 

This restriction is in terms of the reduced interaction and the change to the reduced interaction. 

5V is defined in Eq. (p^, which makes it clear that since V{K') can be expanded in powers of g^,, so 
can 5V: 

oo 

5V = Y.gl5V^'\ (24) 

t=2 

We refer to JV^*^*-' as the ©(f/^, ) change to the reduced interaction, although for convenience the coupling 
is factored out. Note that 5V^^^ is a function of A and A'. 
Now Eq. ( ^3|) can be expanded in powers of 5^ and g^, : 

OQ 00 00 

5: 5l ( A) - ^ 9l ( A') = ^ gl SV'-'y (25) 

t=l t=l t=2 

This equation is a bit tricky to use because it involves the couplings at two different scales. To see how 
they are related, consider the matrix element of Eq. (Esh for gi —* 5233: 



(32.93I V{A) |.gi) - (52331 ^(A') \gi) = (5233I SV . (26) 
According to the definition of the coupling, this equation implies that 



mTT^P+6^^\pi-P2-P3)^92g3\v\gi:ir] '(32331-5^151) . (27) 



Since V^^'^ changes particle number by 1, inspection of Eq. ( p^ reveals that (5233I \gi) is 0(5^,); so 

g.^g.,+o{gl,). (28) 

This implies that 

00 

3a=3v+E5;G^(A,A'), (29) 



where the Cg's are functions of A and A'. For an integer t > 1, Eq. (29) implies that 



3l =.i+E5A?'^*-(A'^'), (30) 



s=2 
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where the Bt^s's are functions of A and A', and can be calculated in terms of the Cs's by raising Eq, 
to the t^^ power. 

We substitute Eq. ( ^o|) into Eq. ( p5|) and demand that it hold order-by-order in g^,. At 
(r > 1), this implies that 

r-1 

V^''\A) - T/('^)(A') = SV^''^ - J2 Br-s,sV^''''\A), (31) 

s=2 

where SV^^"^ = 0, and we define any sum to be zero if its upper limit is less than its lower limit. The 
cutoff-independent parts of V'^^^A) and V^'^\A') cancel on the left-hand side (LHS), leaving 

r-1 

y^2(A) - V^^-^A') = JFf' -J2Br-s,sV^^~^HA). (32) 

s=2 

This equation can be used to derive the desired recursion relations. 

In the remainder of this section, we summarize the results of Appendix B, which contains a derivation 
of the recursion relations that define the IMO. Appendix B is an extension of Appendix D of Ref. [Q to 
the case of pure-glue QCD. 

Recall that momentum conservation implies that any matrix element {F\ V{A) \I) can be written as 
an expansion in unique products of momentum-conserving delta functions. This means that an arbitrary 
matrix element of Eq. (|3^) can be expanded in products of delta functions and thus is equivalent to a set 
of equations, one for each possible product of delta functions. Given approximate transverse locality, each 
of the resulting equations can be expanded in powers of transverse momenta. Matching the coefficients 
of the powers of transverse momenta on either side of these equations allows us to rigorously derive the 
following results (see Appendix B for details). 

First, the cutoff-dependent part of the C(g^) reduced interaction is given in terms of lower-order 
reduced interactions by 



{F\V^'2i^)\I) 



s=2 



(33) 

A terms 



where "A terms" means that the RHS is to be expanded in powers of transverse momenta and only the 
terms in the expansion that are proportional to powers or inverse powers of A contribute. Recall that 
SV^''^ is defined by Eqs. (0) and (H). 

Second, the cutoff-independent part of V^^^A) has two contributions: a four-point interaction with 
no transverse-momentum dependence, and a three-point interaction that is linear in transverse momenta. 
If there are no such contributions to V^^''{A), then it is completely determined by Eq. (^3|). Third, 
the coupling runs at odd orders; i.e. Cg is zero if s is even [see Eq. (p9|)]. Fourth, there is no wave- 
function renormalization at any order in perturbation theory in our approach because this would violate 
the restrictions that we have placed on the IMO. 

The fifth and final result from Appendix B is that the cutoff- independent parts of the 0{g2) and 
0(5^+^) reduced interactions for odd r > 3 are determined by the coupled integral equation^^ 

^'^It is very difficult to prove that integral equations of tliis type have a unique solution; so we simply assume that it is 
true in this case. 
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{F\V^^^ \I) 



1 



1 



r+1 



r+1, 2 



3— point 



{F\ <5F('^+2) |/) - Br+2-s.s {F\ V^('-+2-«)(A) |/) 

s=3 
r+2 

{F\ SV^-+^^ \I) - ^ i?r+3-.,. {F\ F('-+3-^)(A) |/) 



5 = 3 



pj^ term 
4 — point 



p term 



(34) 
(35) 



To use these equations, the right-hand sides have to be expanded in powers of transverse momenta. Only 

(r) 

three-point interactions that are hnear in transverse momenta contribute to V^j , and only four-point 

(r+1) 

interactions that are independent of all transverse momenta contribute to Vqj . 

These equations are coupled integral equations because both Vq^ and Vq-^'^^ appear on the RHS 
of Eq. ( |3^ ) inside integrals in , and Vq^^"^ appears on the RHS of Eq. ( |35| ) inside integrals in 

Sh inside integrals in 5V''^'^^\ but 



^V"('-+3) . It would seem that Vq^^'^^ also appears on the RHS of Eq, 



VJjj cannot couple to ' to produce a transverse-momentum-independent four-point contribution 



to SV^'^'^^K This is because the cutoff function Ti"^'^ ' vanishes when the intermediate state is massless 



and all external transverse momenta are zero. This means that since we specified V^^^ and vif m 



Subsection |2.3.4| , we can use Eqs. (^ and ( |35| ) to solve for Vq^ and V^^' simultaneously, and V^-^' and 



CI 

r(5) 



Vqj' simultaneously, and so on. Note that before we can use these equations to solve for V'^j'' and Vqj 
simultaneously, we must first use Eq. ( |33| ) both to compute V^^^ (A) in terms of lower-order interactions 
and to express Vq-^^\K) in terms of lower-order interactions and ^(''^(A). 



(r+1) 



3 The Basis for the Expansion of Physical States 
3.1 Preliminaries 

In the remainder of this paper, we use the results of our renormalization procedure to compute the 
physical states of pure-glue QCD and their masses. The states will be eigenstates of Al^(A). Since our 
cutoff preserves translational covariance and covariance under rotations about the 3-axis, we would like 
the states to also be simultaneous eigenstates of the generators of these symmetries, but this is impossible 
because translations do not commute with rotations. However, a rotation about the 3-axis separates into 
a part that rotates the centers of mass of states and a part that rotates states' internal degrees of freedom, 
and translations do commute with these internal rotations. 

To be precise, we define jTa to be the generator of rotations about the 3-axis, and to be the part of 
jTa governing gluons' momenta in the center-of-mass frame and spin polarizations. A^^(A), 7^+, V^, and 
are a set of commuting observables; so an eigenstatc of A^^(A) can be labeled by their eigenvalues. 
We choose to write a physical state as [^'■'"(P)), where P is the three-momentum of the state, j is the 
eigenvalue of for the state, and n labels the mass eigenvalue of the state (n = 1 has the smallest mass, 
n = 2 has the second-smallest mass, etc.). Note that because [^'•'"'(P)) will be determined by A^^(A), it 
will implicitly depend on A and . 

An examination of the matrix elements of A^^(A) leads us to believe that the light-front infrared 
divergences (see the discussion below) will not cancel unless the physical states are color singlets, although 
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we do not have a rigorous proof of this. Therefore, we assume that the physical states are color singlets. 
Using the Fock-space expansion for the identity operator (see Appendix A for our light-front pure-glue 
QCD conventions), we can expand a physical state in terms of the number of gluons: 



|vI/^"(P)) = 1 |vl/J"(P)) ^ ^ 1 D,D2{gig2\^^'\P)) 



9192) 



(36) 



where there is no one-gluon component because there is no color-singlet gluon. We neglect contributions 
to the states with more than two gluons, which is a severe approximation. This means that we are 
approximating all physical states as relatively simple single-glueball states, and j is then the projection 
of the glueball's spin onto the 3-axis. From now on, we refer to the approximate physical states simply 
as glueball states. 

Any eigenstate of A^^(A) can be written as a superposition of cross-product states: 



(37) 



where |T„_i) is the color part of the i*"^ contribution to the state and r^"(P)y is the momentum/spin 
part. 

The free state \g1g2) is defined as a two-gluon state in which one particle has quantum numbers pi, 
si, and ci, and the other has quantum numbers p2, S2, and C2. \g1g2) is symmetrized with respect to 
particle labels. Thus we can write 

1 



15152; 



V2 



l5i;52) 



152;. 91) 



(38) 



where the semi-colon in |(?i;.92) indicates that particle #1 has quantum numbers pi, si, and Ci, and 
particle #2 has quantum numbers P2, S2, and 02- Then we can write the wave function as the sum of 
two permutations: 



1 



(gi52|vI'-'"(P)) = ^[(5i;52| + (52;gi|J |vI'-'"(P)) 



(39) 



This implies that the wave function (51(72 |vE'-'"(P)) is symmetric under exchange of particles 1 and 2. 
Using Eq. (|37|), we can write the first permutation as 



(gi;g2|M/^"(P)) = ^ (ci; C2 |T„,,) (pi, si;p2, S2 Tf (P) 



(40) 



(ci;c2|T„^j) and (pi,si;p2,S2 



r^"(P)\ are the color and momentum/spin wave functions for the first 



permutation of the i contribution to the total wave function. 



3.2 Calculation of the Color Wave Function 

In our two-gluon approximation, the fact that the glueball states are color singlets uniquely determines 
their color wave functions. Using the completeness relation for color states, the color part of the i^^ 
contribution to the glueball state can be written 

l^n.i) = ^ |ci;c2) (ci;c2 |T„,i) . (41) 

C1C2 
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Let Y be the unitary operator that rotates two-gluon states in color space. Then under a color rotation, 
we have 

|T„,j) ^ |T^_,) = ^F|ci;c2) (ci;c2 |T„^,) 
C1C2 

= $]|c'i;4)(ci;c2|T„,,) , (42) 
C1C2 

where the primes on the basis vectors denote rotated basis states. Since |T„^i) is a color singlet, it does 
not change under a color rotation: 

\^n,^) = = E l<;^2> (C1;C2 |T„,.) . (43) 

C1C2 

A gluon color-rotated state is given by 

Ic'i) = U |ci) = J2 |C2) (C2I U |ci) = ^==-1 1^2) , (44) 

C2 C2 

where 

[;^g-»e.i#=i (45) 

and 

^C2^C3 = -*r^'^'^- (46) 

Here 9c is a vector of — 1 real numbers that parameterize the color rotation. Since 

|ci;4) = |c'i)0|4), (47) 



Eqs. (Ua) and (44) imply that 



As a guess, let us try 



^ (ci;c2 |T„,j) C/c3CiC/c4C2 |c3;c4) • (48) 



ci;c2 |T„,,) = 4^ci,c2, (49) 



N 

where N is some constant. Then Eq. (|4^) becomes 

= ^ E t^c3c,J7c,cJc3;c4). (50) 

C1C3 C4 

Note that 

([/^l)cc' - {U^)cc' = UX = f/c'c (51) 

Then 

\'^n,i) = E '^'^S.Ci |C3;C4) 

C3C4 



(ci;c2 |T„,i) |ci;c2) 



C1C2 



. (52) 
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Since this is true, our guess for a color-singlet wave function was correct. This solution is unique psl; so 



we can drop the subscript i in Eq. ([40[). N is determined by normalization to be y/ — 1. 
Now Eq. ( ^9|) becomes 



(9192 = , ^ 5c, c. 



(53) 



Since the IMO conserves momentum, the factor in brackets must be proportional to a momentum- 
conserving delta function. Thus we can write 

(gi52|*^"(P)) - 2!(16^3)i^(3)(p_p^ (54) 

V2(^c - 1) 

where $^"^2 (f^i'Ps) is the momentum/spin wave function with the momentum-conserving delta func- 
tion removed, and all the extra factors in this equation are present to simplify the normalization of 
^^"^^(Pi'Pa)- We must solve for ^^"g^ (pi,p2), and this equation indicates that it is symmetric under 
exchange of particles 1 and 2. 

3.3 Jacobi Variables 

Using Eq. (|5^), the Fock-state expansion for a glueball state in Eq. (|3^ ) becomes 

Vl6^3 ^2(7V2 - 1) ^^^^^^^^ J 

X *ir.,(Pi,P2)|3i52), (55) 

where momentum conservation implies that P2 = P ~ Pi- It is useful to separate the motion of the center 
of mass of the state from the internal motions of the gluons. To do this, we change variables from pi to 
the Jacobi variables x and k±: 

pi = (xP+,xP^+k^), 

P2 - {[\-x\P+,[\-x\P^~ki_). (56) 

Here x is the fraction of the total longitudinal momentum that is carried by particle 1, and fc^ is the 
transverse momentum of particle 1 in the center-of-mass frame. We only display the longitudinal and 
transverse components of the momenta. (Since the glueball state is a superposition of free-particle states, 
and since the momentum of a free gluon satisfies p^ = 0, the minus components of the momenta of gi 
and g2 are constrained to be given by p~ = Pi±/pt ■) 
In terms of the Jacobi variables, the glueball state is 

= /I 3 Mm U ^ Sc,^cJ^f^eix-e)eil-x-e) 
Vl67r3 ^2(iV2 - 1) J ^x{l - x) 

X $f;,Ja;,fcj.,P) g{x,k±,P;si,ci)g{l-x,-k±,P;s2,C2)) , (57) 

where we explicitly show the dependence of the RHS ket on the Jacobi variables and the total momentum. 
We can use boost covariance to show that $^"^^(0;, k±, P) is independent of P. To do this, we note that 
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under a longitudinal boost, the longitudinal momentum of each particle in pure-glue QCD (whether the 
particle is point-like or composite) transforms according to 

p+ e>+, (58) 

where v is a. boost parameter p^ , p5| . Under a transverse boost, each particle's transverse momentum 
transforms according to 

P±^P±+P^v±, (59) 

where v±_ is a boost parameter. This means that if we apply a boost operator to both sides of Eq. (p7|), 
we find that 

1*"^^')^ - /I 3 Mm ^^ ^ 5.,. / ^!^^(x - e)^(l - x - e) 



g{x, k^,P'; si,ci)g{l - x, -k^,P'; S2, 02)^ , (60) 



where the boost takes the glueball's momentum from P to P' . Note that the boost does not affect the 
wave function, only the kets. Since Eq. ( p7| ) holds for all P, it holds in particular for P': 



X ^l^,^{x,k^,P') g{x,k^,P';si,ci)g{l-x,-k^,P';s2,C2)). (61) 

Eqs. (|60| ) and (61) contradict each other unless the wave function ^l^s^ix, k±, P) is independent of P: 

$f;,^(a:,fc^,P) = $f;^,(x,fc^). (62) 

Thus we can write |4'^"(P))) as 

/I 3 /0!M2 U E K,cJd^k^dx I ^e$i^,jx,fcj|gig2), (63) 

V167r3 V2(iV2 - 1) ^^^^^^^^ J ^a:(l - x) 

where 9^ = 9{x — t)0{\ — x — e). 

3.4 The Momentum and Spin Wave-Function Bases 

To solve for ^^)' expand it in a complete orthonormal basis for each degree of freedom. Since 

^{"^g^ix, k±) is symmetric under exchange of particles 1 and 2, 

n^s, i^, k±) = (1 - 2:, ~k^). (64) 

To use this, we define 

4 

<^Z,{x,h) = 5]x?^^^f^r(^'^^)' (65) 
9=1 

where the spin wave functions are 
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Xs''" = [^s,,iSs^,i + Ss,,iSs2,i] , 

xl'"^ = -^[Ssi,iSs^^i- Ss,,iSs2,i], (66) 
where s — —s, and the momentum wave functions satisfy 

ni'\x, fc_i) = f7f (1 - X, -fc_i), 

nl"{x,k^) = r!f (l-a;,-fci), 

ni'\x,k^) = -r!f (l-a;,-fci). (67) 

Note that 

E^'^^'^--^,,,. (68) 

SlS2 

We define k = \k±\, and we define the angle (p by 

fcj^ = fc cos a; + fc sin (/) y. (69) 

We expand the momentum wave function in complete orthonormal bases: 

oo oo oo 

n^{x,h) = Y.ll E RluaL\'\^m^'\k)AM), (70) 

1=0 t=0 a=-oo 

where ^[''^(a;), Ti'^\k), and Aa(0) are basis functions for the longitudinal, transverse-magnitude, and 
transverse-angular degrees of freedom, e and d are parameters that govern the widths of the longitudinal 
and transverse-magnitude basis functions, respectively. We can adjust these widths to optimize the bases. 
Note that if we do not truncate the sums in Eq. ([70|), then the R'JitaS depend on e and d such that 
51^" (a;, kj_) is independent of e and d, although we do not indicate this dependence explicitly. 
We define the transverse-magnitude basis functions T^ '^\k) by 

t 

ri''Hft) = de-'=''^'E'^*.«^'^'' (71) 

where the crt,s's are constants that we have computed numerically using the Gram-Schmidt orthogonal- 
ization procedure and are such that 

dk k T^/^ {k)T^f (k) = St,t' ■ (72) 



Adjusting d allows us to adjust the width of the Gaussian weight function in Eq. ([TlJ). d cannot be 
allowed to pass through zero, because when d is zero, T^'^^k) is zero. We choose d to be positive without 
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loss of generality. When doing numerical computations, we work with a dimensionless form of these basis 
functions: 

T,{kd)^^T,^''\k). (73) 



Under exchange of the two particles, k is unaffected; so Tf'^\k) is unaffected. 

'I 



We define the longitudinal basis functions ll'f' {x) by 



L|'^)(.x) = Kl-x)r^A^)„x™, (74) 



m=0 

where 



^(e) ,y-„. 1 r(l+4e + ; + m) //!(l + 4e + 20 

'■"^ ^ ' m\{l-m)\ r(l + 2e + m) y r(l + 4e + ' 

These definitions imply that 

dxLf{x)Lf{x)^8i,v, (76) 

as long as e > —1/2. (If e < —1/2, then the state is not normalizable.) We can adjust the width of the 
weighting function in Eq. (^3) by adjusting e. It is often more convenient to work with 



-L^:^[x)^^^=Lf\x). (77) 
V 2:(1 — 



Under exchange of the two particles, x ^ \ — x^ which means that L^f^ (x) (— l)'i|'^'' (a;). 
We define the transverse-angular basis functions Aa((/)) by 

A,{^) - ^e-^ (78) 



/27r 
and then 

27r 

d^Al,{i,)Aa{<i>)^da,a'- (79) 

These basis functions are useful because they are eigenfunctions of £3 , the part of the generator of 
rotations about the 3-axis that governs gluons' momenta in the center-of-mass frame. Under exchange 
of the two particles, —>(/) + tt, which means that Aa(0) (— l)"Aa(0). 
Using the above definitions of the bases, |^-'"(P)) can be written 

|*^"(P)) = ^/g^3 E ^^uc.KnaX'r' J d^kML^t\^)Ti''\k)Aa{^) \g1g2) ■ (80) 



sxS2CiC2qlta 



Since the glueball state is symmetric under exchange of the two gluons, the behaviors of the spin and 
momentum wave functions under exchange of the two gluons imply that if g = 4, then I + a must be odd, 
and if g 7^ 4, then I + a must be even (so that the spin and momentum wave functions have the same 
symmetry under exchange). 
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3.5 Rotations about the 3- Axis 



We want to ensure that [^'■'"(P)) is an eigenstate of with eigenvalue j. The action of jTs on the state 
is given by 

1 1 



siS2CiC2qlta 

d 



dp 



W2^ 



1± 



dpi 



2± 



dp. 



+ Sl + S2 



2± 



Lf\x)Ti''\k)A,{<i>) 



X 15152) 



(81) 



Using the definitions of the Jacobi variables, we can separate the center-of-mass and internal degrees of 
freedom: 

1 1 



J3|*^"(P)) 



\/l67r -yJli^N^ ry g^52CiC2Q/^a 



R. 



-'ci,C2-"'g;t, 



Sl+ S2 



d^k±dxOe 



Ll'''ix)T^'\k)A4cp)}\g,g2). (82) 



This implies that the action of on the glueball state is given by 



J3^|*^"(P)) 



Vi6^ vmi - 1) 



. d 
dcp 

1 1 



S2CiC2qlta 

L\'\x)T,^''\k)A4cl>)}\gi92) 



V V C / 5l52ClC2Qtta. 

X T,^''\k)A,{<P)\g,g2). 



itaAq j 



d^kj_dx9e [a + si+ S2] L^f' {x) 



(83) 



Since [^'■'"(P)) is an eigenstate of with eigenvalue it must be the case that a + si + S2 = j for all 
values of a, si, and S2 that contribute to the sums in this equation. This implies that for some set of 



coefficients P^J^, 



This means that |^-'"(P)) can be written 



|*i"(p)) = ^P^^pjg,Z,t,j), 
qlt 



(84) 



(85) 



where 
|9,^,i,i) 



— E 



d'k^dxe,5,,,,,x\'''Lf\x)T^'^\k)Aj_s,-sM \9192) . (86) 



We are going to calculate the matrix elements of the IMO in the \q. I. t.j) basis and diagonalize it. 
Since M'^{K) commutes with we can do this for each value of j separately. The diagonalization 
procedure will yield mass eigenvalues and the P^"^ coefficients. As long as the coefficients satisfy 



A pi'"'* pi" _ X A 

"id' ^qlt ^qlt ~ '^j,3'"n, 
qlt 



n' ) 



(87) 
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the glueball state will have a plane-wave normalization: 

($/"'(P')|5-J»(P)) = 167r3p+(5(3)(p _ P')Sjj,6n,n'- (88) 

Due to the symmetry of the glueball state under the exchange of particles 1 and 2, the basis state 
\q,l,t,j) is zero if I + j is odd and q 7^ 4, or if ^ + j is even and q ~ 4. To take advantage of this, we 
consider only the subspace in which I + j is even if g ^ 4, and 1 + j is odd if g = 4. 

3.6 The Eigenvalue Equation 

The IMO's eigenvalue equation is 

M^{A) |\E'^"(P)) = M2 |*J"(P)) , (89) 

where the lower-case subscript on indicates that it is an eigenvalue of A^^(A) rather than an eigenvalue 
of A^free- (An example of an eigenvalue of A4%^^ is M^.) 

In the basis that we have defined, the eigenvalue equation takes the form 

M\A) 5] R^^ \q, I, t,j) = Ml J2 %l l9. ^. td) ■ (90) 

qlt qlt 

If we project {q', I', t',j\ onto the left of this equation and use the identity 

{q',l',t',j\q,l,t,j) = 16n^P+6^^\P-P')Sg,g,6i,i'6t,t', (91) 
then we find that the eigenvalue equation in our basis is 

^ {q',l',t',j\M'^{A)\q,l,t,j) ^jn _ ..2r>jn .q„x 

^ 167r3P+5(3)(P-P') ~ ^^^n^g'i'f ly^J 

qlt ^ ' 

4 Calculation of the Free-State Matrix Element of the Invari- 
ant-Mass Operator 

4.1 Preliminaries 

To solve the eigenvalue equation, we must calculate the matrix elements {q' ,1' ,t' ,j\M.'^{h-)\q,l,t,j). 
These can be written in terms of the free-state matrix element (3^52! ^^^(A) \gig2)-, which is specified 
by our renormalization procedure. We are going to calculate this matrix element to second order in 
perturbation theory. 

Before continuing, we would like to point out three simplifications that we repeatedly use in this 
section. To make these simplifications clear, we note that the matrix elements (q', i', i\ A^^(A) |g, t, j) 
can be written as integrals of wave functions times the free-state matrix element {g'ig2\ {h) \gig2)- 
To make the first simplification, we observe that {g'ig2\ M/^ {A) \g1g2) has step functions on the particles' 
longitudinal momenta, and these step functions are redundant because they also appear in the integrals 
that define {q' ,1' ,t' ,j\M.'^{K) \q,l,t,j). We drop these step functions in the formulas that we present 
for {g'^g'2\M.'^{K)\gig2). To make the second simplification, we point out that the matrix elements 
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{q\l\t' , j\ M.'^{K) \q,l,t,j) are symmetric mrder exchange of the two initial-state gluons and also under 
exchange of the two final-state gluons. Because of this, when we are computing {g'ig2\ A^'^(A) 13152), we 
combine terms that are related by exchange of the two initial-state or two final-state gluons. Finally, due 
to the color-singlet nature of the glueball states, certain parts of {g[g2\M'^{A) \g1g2) do not contribute 
to {q',l',t',j\M^{A)\q,l,t,j), and we drop these terms. 

We begin by defining Jacobi variables for the final state of the free matrix element: 

p'2 = ([1 - x']P+, [1 - x']Px - fcl). (93) 

Then using the definition of the IMO in terms of the reduced interaction in Eqs. (P, (|^), and (p^, 

{gW2\MHA)\g,g2) = (^iffil A^Lc IffiSs) + (fflssl ■^f„t(A) |5iff2) 

= m! {gW2 1 5152) + fff e-^"'^" {gW2 1 Igm) 

+ 5>-^"^"(ff'i52l^(^'i](A)|5i52), (94) 

where the free masses of the final and initial states are given by 

U2 

Ml 



x'{i~x'y 

Mf = — (95) 
x{l — X) 

4.2 The Cutoff-Independent Part of the Reduced Interaction 

Based on the definition of V^^^^ in Eq. (p^, we find that 

4 

(3I32I 15132) = 167r3p+5(3)(P _ P') ^ ^, j:g[g'Mglg2:^^. (96) 

1=1 

We can divide this matrix element into contact (momentum-independent) and instantaneous-exchange 
interactions: 



where 



(5152! ^ci 15132) = (g'i52l ^ci l5i32)coN + (ffi52l ^ci l3i32)iN , (97) 

{9W2\vg^ \9i92)coN = 167r3F+5(3)(F-F')^l!ii',2' <^5i32kl3i52^i 

= 327r3p+j(3)(P- P')r^^'^7^^^^^ - K'.rs'.h.M) , (98) 



and 



k9x92\ ^(^?^ l3i52)iN = 167r3p+<5(3)(F - P') ^ ^ ^g\g'2\v\gYg2-\ i 



^(2)|„. _ ifl„3p+x(3)/p 

i=2 

327r3p+,5(3)(P - P')B{\x - x'\ - e)r^"'i"r^"^"(5,,,,;5. 



^ -j^ix + x')(l - X + 1 - x'). (99) 



(a; — x') 

(The remaining term in Eq. ( p6|) vanishes because the |g, Z,t,j) states are color singlets.) The contact 
interaction is displayed in Figure 1 and the instantaneous-exchange interaction is displayed in Figure 2. 
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1' 

2' 

(2) 

Figure 1: A diagrammatic representation of {g[g2\ V^j |3i52)coNi two-gluon contact interaction. The 
numbers label the particles. 

1 JLOJJLOJLOJJl 1' 




2' 

(2) 

Figure 2: A diagrammatic representation of (<7i52l^ci |5i52)in; the two-gluon instantaneous-exchange 
interaction. The numbers label the particles. 

4.3 The CutofF-Dependent Part of the Reduced Interaction 
4.3.1 Preliminaries 

According to our recursion relation for the cutoff-dependent part of the reduced interaction, Eq. (p3|), 
{91921 ycuW \9192) is given by 

{9W2\vgliA)\g^g,) = {9W2\SV'-'^ \g^g,)\^^^^^^. (100) 
Using the definition of the 0(5^,) change to the reduced interaction in Eqs. ( p^ and (|2^), we find that 
MI^(^?](A)|5i<?2) = i5](ff;5^|FW|if)(i^|T/«|5i52)Tf'^'^(^^,i^,/)|^^_. (101) 

K 

The intermediate state can be either a one-particle or three-particle state. The contribution to the 
eigenvalue equation from the one-particle-intermediate-state part is zero because the \q,l,t,j) states are 
color singlets. This means that 



(5'i52l ^cd(A) \g1g2) - ^/ DsD^D, {gW,\ V^'^ I533455) (ff35455| V^' 



^ '9152) 



X Tt^''\F,K,I)\^^_, (102) 

where Mj^ is the mass of the state 1533455). Substituting the definition of y^^^ in Eq. ([ij) into this 
equation, and simplifying, we find that 

(5I52I vgliA) 15152) = (5i52l vgliA) |5i52)se + (5i52l vgl{A) |5i52)ex : (103) 

where the self-energy interaction is given by 

(5i52l^^CD(A) |5i52)sE = P+^IQtt')' I D^D^D.T^^^^'^F, K, I) 62,, S^' ,5 PiMf3*v3,\ ^^^^^^ , (104) 

the exchange interaction is given by 
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1' 



2 Trrr^TTTTTTJ 2' 

(2) 

Figure 3: A diagrammatic representation of (5'iff2l ^cd (^) |5'i52)se' self-energy interaction. The 
numbers label the particles. 



1' 



2' 



Fi gure 4: A diagrammatic representation of (3152! ^^d(''^) |5i52)ex' exchange interaction. The num- 
bers label the particles. 




(5'i52|l^CD(A)|5i.92>Ex-2P+2(167r3)2y D^D^D^4''^'''\f, KJ) 5^,^ -^i'^^ /3i34/32* 45 L , (105) 



and we have defined 



and 



13,^^ = e{pt)e{p+)e{pl)5^^\p,-p, -pk) UJ9k\v\g^-}r 



(106) 



= I&n^pt5'^^\p,-Pk)5s,,sjc,,c,- (107) 
The self-energy interaction is displayed in Figure 3 and the exchange interaction is displayed in Figure 4. 

4.3.2 The Self-Energy Interaction 

After some simplification, the self-energy interaction takes the form 

cPpz^dp'l 



(3l52lK^?](A)|ffi52; 



SE 



1 



-^e{p+ - eP+)^9{p+ - eP+) 
P3 



Pi 



A terms 



(108) 



where P4 = Pi — Ps- At this point, it is useful to change variables from p^ to Jacobi variables y and f±: 

P3 = {ypt,ypi± + r±). (109) 
Then the free mass of the intermediate state is 
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+ 



(110) 



x{l - x) xy{l - y) ' 

and after a lot of additional simplification, 

{9W2\ V^^k^) |5iff2>sE = ^\/| [A' - A'^] <5i,i, ,52,2' 9{x - 2e) 

^■"^ V ^ L -^^J A terms 

Recalling that "A terms" means that we are to expand the RHS of this equation in powers of transverse 
momenta and keep only the terms that are proportional to powers or inverse powers of A, we find that 



log a; - log e - 



(111) 



{gW2\ ^(^S(A) \gi92)sE = £v|a2 ^1,1, 62,2' e{x - 2e) 



logx-loge- 



(112) 



This result shows that our cutoff violates cluster decomposition. The evidence is that the self-energy 
depends on. x = p'^ / {p^ + J^J)) even though p2 is the momentum of a spectator (see Figure 3). 



4.3.3 The Exchange Interaction 

For the exchange interaction, the free mass of the intermediate state is 

x' x — x' 1 — X 

Then the changes in free mass are 

fc'2[l - a;]2 + _ xf - 2kk'{l - x'){l - x) C0S7 



where j = (f) ~ cj)'. 
Using the identity 



{1 - :r:'){x ~ x'){l - x) 
fc^x'^ + fc'^x^ - 2kk'xx' cos 7 
xx'{x — x') ' 



e^s-ki. = -^kse''^, 
v2 



the exchange interaction becomes 

{9W2\ ^(^S(A) |5i52)eX = 647r3p+5(3)(P - p')fc.<fcc.c'. 



e{x -x' -e) 



X — X' 



-2A' '^AfkAik _ g-2A "AfkA/k 



i,m=l 



A terms 



where 



:^sis[Ss,y, j-fce'"^'^ - k'e''"'''] \k'e-'''^'^' \1 - x] - ke-'''^'^[l - x']] 

■ X \ X J V J 



(113) 



(114) 



(115) 



(116) 
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^' ^ { X ) { 1 — X 

Q^'''^ = 25s,y^S,^y^j-^-^ {x{l - x)k'^ + x'{l - x')e - kk' [x'{l - x) + x{l ~ x')] COS7} 



1-x' 
1 - 
~x 



, 1 — r' 



g(2,2) ^ J ^, S2s' (x/c'e'"^'^' - x'ke''"f'] Ike-'''^'^ 
q(2,3) ^ 324(5.1,.;^^:^ [xk'e-'''^^' - x'ke~'''^^^ jfc'e 

' ■ L X J 1 — a;' 



g(3.3) = s2s;<5,,.4{^fcV*'i^'-fce-*^i^}<j/fcV^^^'-i-^fce'^^4. (117) 



In Eq. (116), if we expand everything multiplying the delta function in powers of transverse momenta, 
the lowest-order terms from the two exponentials cancel and leave two types of terms: those proportional 
to inverse powers of A and those proportional to inverse powers of A'. We can isolate the terms that are 
proportional to inverse powers of A without altering the cancellation of the lowest term by replacing the 
first exponential with a 1: 

3 



^FK ^IK 



i.m—1 



4.4 Combining the Interactions 

In order to get the infrared (e 0) divergences to cancel, it is useful to combine the interactions in a 
particular manner. From Eqs. (p3), (p7|), and (|103|), 



{9W2\ ■M'(A) 13132) - {gW2\M^A) Igm)^^ + (siffil M\A) |gi52)coN + {9W2\M^A) 1.91,92)1 



)sE + (.9^.9^1 >''(A)|gig2)Ex. 



{gW2\M'{A) |.9ig2)sE + (siffil M^{A) \gig2)Ex , (US) 



where 



{gW2\M\A)\gig2)KE = Mf {g[g'2\9ig2) , 

{g[g'2\ M\A) |ffi52)coN ^ .9>-^"^" (ff'iSsI ^c^i^ l.9iff2)coN > 

{gW2\M\A)\g^g2\^ = .9>-''"^" (5'i52l ^c^?' I.9i52)in . 

(31521 M'(A)|5i52)sE = 5>^^"^"(5'i52l^c^S(A)|gi.92)sE, 

(gl52l>''(A)|.9i.92>EX = S^"^"^" (ff'i52l ^(^'i](A) |.9i.92>ex ■ (120) 



We break the instantaneous interaction into two parts: a part that is "above" the cutoff, i.e. a part 
that would vanish if we took A — )■ 00, and a part that is "below" the cutoff, i.e. a part that would survive 
if we took A — > 00: 
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{g[g'^\ M'{A) \g1g2U = (ff'i^sl M\A) Igm)^^ + (ff'iSsI -^'(A) \gig2)f^ , 



(121) 



where 



/IN 



(gl52|A^'(A)|.9ig2)fN = (l-e-^^"^--^^-)(g'i.9^|A<^(A)l5i52> 
(g'i52l-M'(A)|5i52)fN = e-2^"'^--^-(ffl52l-M'(A)|5iff2)iN- (122) 
Next, we break the self-energy and exchange interactions into finite and divergent parts: 

{g[g',\ M\K) |gi.92)sE = (siffsl M^K) |gi52)sE + {g'MM^k) |5i52)se . 

{g'M M^A) |.9ig2)Ex = (s'lffal -^'(A) |5i32)ex + (siffil -^'(A) |.9i.92)ex > (123) 



where the divergent part of the self-energy interaction consists solely of the term containing the loge, and 
the divergent part of the exchange interaction consists solely of the term containing Q^^'^\ Finally, we 
define an interaction that is a combination of the instantaneous interaction "above" the cutoff and the 
divergent part of the exchange interaction: 

(gW^lM^A) |5i52)iN+EX = (siffil ■^'(A) \gm)t^ + (ff'iffil A^'(A) |gi52)Ex • (124) 

Then 

{g[g',\ M\A) \g,g2) = {g'M M\A) \gm)KE + (ff'iffil ■^^(A) 131.92)^^ + {gW2\M\A) |5iff2)coN 
+ {g[g'^\ M\A) \gig2)lx + (g'lfhl M'iA) \g,g2) 

IN+EX 

+ {gW2\M\A)\g,g2) 

IN 

+ {gW2\M\A)\g,g2)^^. (125) 

Perry showed that with a suitable definition of long-range interactions, a renormalization method that 
is similar to ours yields a logarithmically confining potential for quark-antiquark bound states at 0{g'^) 
[0 . His calculation uses sharp step- function cutoffs and is based on an analysis of the part of the two-body 
interaction that is most singular in the limit in which the exchanged gluon has infinitesimal longitudinal 
momentum. The corresponding part of our interaction, which is contained in ((7^.92! M'^{A) |5i(72)in+ex 
and (5132I -^'^(A) |3i52)fN' similar to what Perry found in the quark-antiquark case. However, to 
determine whether or not our interaction is truly confining, we would have to do a careful analysis of the 
complete two-body potential, not just the most singular part. This analysis would be complicated by the 
smooth cutoff that we employ, and we leave it for future consideration. 



5 Calculation of the Matrix Elements of the Invariant-Mass Op- 
erator in the Basis for Physical States 

The matrix elements (9', l\t' , j \ A^^(A) \q,l, t,j), which appear in the eigenvalue equation, can be divided 
into contributions corresponding to the different terms in Eq. (|125|): 
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{q', I', t',j\ M^A) \q, I, t,j) - {q', I', t\j\ M^{A) \q, I, t,j)^^ + {q' , I', t',3\M\k) \q, I, t, j) 



/SE 
■\F 



+ {q', I', t', j| M\K) |g, /, j)coN + J\ M\K) |g, /, j) 



EX 



+ {q',l\t\j\M\K)\q^,t,3) 
+ {q',l',t',j\M'{A)\q,l,t,j) 



IN+EX 
D 

SE ■ 



+ {q',l',t',j\M'{A)\q,l,t,j)f, 



(126) 



In this section we express these different contributions in terms of integrals that can be computed numer- 
ically. These integrals fall into two classes: two-dimensional and five-dimensional. We treat each class of 
integral separately. 



Each of the terms in Eq. (126) is proportional to the plane-wave normalization factor 16n^ P'^ S'^^^ [P~ 
P'). To make the remaining equations that we present simpler, we suppress this factor. We also take 
e ^ in any contribution to {q' ,1' ,t' , \q,l,t,j) that is finite in this limit. 



5.1 The Two-Dimensional Integrals 
5.1.1 The Kinetic Energy 

Using the definition of our basis in Eq. (p6|) and the expression for the free-state matrix element of the 



kinetic energy in Eq. (120), we find that 

<2 



1 

d2 



dxL'jf\x)Ll''\x) I drr^Tt'{r)Tt{r). (127) 



((?',r,t',j|X^(A)|g,/,i,j)KE = V 

" Jo Jo 

Note that the kinetic energy is infinite unless e > 0, whereas normalizability requires only e > —1/2. 

These integrals can be computed with standard numerical integration routines, but the results can 
have large errors for large values of the function indices I, I', t, and t' , due to the oscillating nature of 
the basis functions. It is better to rewrite the integrals as sums that can be computed numerically with 
Mathematica |2^ to any desired precision. Using the definitions of the basis functions. 



{q',l\t\j\M'iA)\q,l,t,j)^ 

t t' 
5Z ^ o't',s'2~ 



r(2e) 



E,(e) ,(e) T{2e + m + m') 

i,ra 2^ '''™'r(4e-Hm-Hm') 



m=0 



m'— 



"^r I 2 



s=0 



(128) 



5.1.2 The Finite Part of the Self-Energy Interaction 



The self-energy interaction conserves each particle's momentum, thus Mf = 7\f|, for this contribution. 



This means that the Gaussian cutoff factor in Eq. ( 120 ) has no effect on the self-energy. For the finite 
part of the self-energy, we use the same method for evaluating integrals that we use for the kinetic energy. 
This yields 



{q',l',t',3\M^{K)\q,l,t,j)l^ 
11" 



dxL[f\x)L\''>ix)xil- x) 



log a; 



12 



(129) 



29 



-Sq,q'St,t' 



47r2 



* ^ m=0 m'=0 



(e) r(l + 2e + m + m') 
r(2 + 4e + rn + m') 



X + 2e + m + m') - %l){2 + 4e + m + m')] ^ , 

where the digamma function ipi^) is given by 



(130) 



dT{z) 
dz 

nz) 



(131) 



5.2 The Five-Dimensional Integrals 
5.2.1 The Contact Interaction 

Using the definition of our basis in Eq. (|6|) and the expression for the free-state matrix element of the 
contact interaction in Eqs. (^) and ( |l20| ), we find that 



{q\l\t',j\M'{A)\q,l,t,j) 



X eik)0{k')Tlf\k')Ti'^\k) I dxdx'e{x)e{l - x)e{x')9{l - x')L\f>{x')L\'''{x)e 



CON 



87r2 



I 



dkdk'kk' 



(132) 



The reader may have noticed that tliis integral is not five-dimcnsional as we have implied, but rather four- 
dimensional. However, when we numerically compute the integrals that have more than two dimensions, 
it is most efficient if we combine them into one integral; so we want their integration variables and their 
ranges of integration to be identical. Thus we increase the number of dimensions of this integral by one 
by introducing an extra integral over 7 ~ — using the identity 

1 



1 



27r 



d76'(7)6l(27r-7). 



(133) 



This will help us to combine this integral with others that contain integrals over 7 that cannot be done 
analytically. 

Since the integration domain of the exchange interaction is restricted so that x > x' when e ^ [note 
the step function in Eq. (118)], we would like to enforce this restriction in the other contributions. Using 
the identity 

l = e{x- x') + e{x' ~ x), (134) 
we break the longitudinal integral in Eq. ( |132[ ) into two parts: 

dxdx'e{x)e{i - x)e{x')e{i - x')L\f>{x')L[''\x)e-^'^'^^^'[e{x - x') + e{x' - x)]. (135) 

In the second term, we let a; ^ 1 — x and — > 1 — x', and then the longitudinal integral becomes 



dxdx' e{x)eil - x)e{x')Q(\ - x')Lf{x')L^f^{x)e 



'^Q{x-x') l + (-l) 



(136) 



Recall that we are restricting ourselves to the subspace of the states |g, in which j is even if 

5^4, and I + ] is odd if q = 4. Then since q = q' for the contact interaction, I and /' must both be even 
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or both be odd for this interaction. This means that the two terms in Eq. (|136D are equaL Thus we can 
write the contact interaction contribution as foUows: 

Nc9' 



{q'J',t',j\M\A)\q,l,t,j) 



CON 



87r3 



where 



D = dx dx' dk dk' d-y k k' e{x) 6(1 - x) e{x') e{x - x') e{k) e{k') 61(7) 6l(27r - 7), 



(137) 



(138) 



and 



C = L^\x')L'f\x)T^f\k')T^/\k) 



(139) 



5.2.2 The Finite Part of the Exchange Interaction 

In order to simphfy the contribution to (g', Z', t', j| A^^(A) \q,l,t,j) from the finite part of the exchange 
interaction, we wish to change variables from to 7 = </> — 0': 

d^k^d^k'^ = dkdk'd(t>d(t>'kk'9{k)9{k')e{<j))9{2Tr-(f>)0{(l)')9{2Tr^<j)') 

= dkdk'djd(j)'kk'9{k)9{k')9{-f + (t)')9{2TT--f-(l)')9{(t)')9{2n -((>'). (140) 

Ah the contributions to {q' ,1' ,t' , j\ A4'^{A) \q,l,t,j) depend on 7 only through dependence on COS7 and 
sin 7. This means that we can use the identity 



to write 



27r-0' p2Tr 

d'-ff {cos '-f, sin j) = / (i7/(cos7, sin7) 



d^k^d^k'^ ^ dk dk' d-y d(j)' k k' 9{k) 9{k') 61(7) 9{2ti - 7) 9{(j)') 9{2ti - (j)'). 



(141) 



(142) 



Inspection of Eq. (116) implies that the integrals in (g', Z', t', j| A^^(A) \q,l, i, j)Ex depend on complex 
exponentials of </>' and 7. However, using the identity 

(i7/(cos7) sina7 = 0, (143) 

where a is an integer, it is possible to trivially do the 0' integral in {q' ,1' ,t' ,j\M.'^{K) \q,l,t, and 
write the remainder as a real quantity with integrals that depend on 7 only through cos 7 and sin 7. The 
result is 



{q\l',t\j\M^{K)\q,l,t,j) 



EX 



1 



1 



1 



A 



FK 



A 



IK 



X 1-e 



-2A"*Ai?A'Ar 



kk' 



where some of the S^y„,'s and si^l's are given by 

q,q q.q o J 



a;(l-x)x'(l-a;') 



(144) 
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and 



S[]l = cos([j-2]7), 

5g = -cos(b--l]7)Kl-a;')+x'(l-x)], 

S[% = ±cos{[j-l]j){x[l-x']+x'[l-x]){x'^ + [l-x'f), 

Si]l = -Lcos07)(l-2x'), 

5g = ^cosi\j-lh)il-2x'){x[l-x']+x'[l-x]), 

Si]l = ^cos([j--2h)[x2 + (l-.x)2], 

Si]l = cos{n)[x^ + {l-xf-2x'{l-x')], 

S^^l = -cos(j7)cos7(a;[l-a;']+a;'[l-a;])(l-2a;[l-a;]-2a;'[l-a;']), 

c(i) _ n 

^fi = sin7sin(j7)(2a;^-2a;^[l + a;']+a;'[l-2a;']+a;[l-2a;' + 4a;'^]), 

Si]l = -Lcos([j--2]7)(l-2ar), 

c(i) _ n 

S[]l = cos(j7)(l - 2a; - 2x' + Axx'), 

S'f] = -cos7COs(j7)(l-2a;)(l-2a;')(a; + a;'-2a;a;'), (145) 



Si'l=SiZ\ , (146) 



The rest of the S^^^, 's and S^^^, 's are given by 



where i = 1,3. The <S'^^g/'s are given by 



c(») _ _ n 

•^1,2 — '^2,1 — 

c(i) _ c(») I 

c(*) _ qii) I 

'-'2,3 ~ '-'1,3 . 

qii) _ aii) I 

'-'2,4 — '^1,4 . .1 

q{i) _ qii) I 

'^3,2 — '-'3,1 . .' 

Sfl = -Sm. (147) 
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5.2.3 The Instantaneous and Exchange Interactions Combination 

Using similar metliods for the contribution to (g', Z', t', j| A^^(A) \q,l,t,j) from tlic combination of the 
instantaneous interaction above the cutoff and the divergent part of the exchange interaction, we find 
that 



9; ^' j)lN+EX 

2 

{x + x'){l - X + 1 ~ x') + 



87r3 
f 



1 
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x{l-x)k''^ +x'{l-x')k'^ 



kk' [x'{l - x) + x{l - x')] cos 7^ 



(i49) 



where 



Wq,q> = Sq^lSq'S COs{-f[j - 2]) + 6q^2Sq' ,2 COs{j[j + 2]) + (5,^3^:3 C0s(7j) + 5q^45q'A COs{'yj). (150) 

Note that the divergences from the two interactions that comprise {q' , I', t' , j\ A^^(A) \q, I, t, j)in+ex 
eel, allowing us to take e — > in this contribution. 



5.2.4 The Instantaneous Interaction Below the Cutoff 

The contribution to {q' ,1' ,t' , j\ M'^{A) \q,l,t, j) from the instantaneous interaction below the cutoff is 
divergent. In this subsection, we extract the divergence, show that it cancels the divergent part of the 
self-energy, and compute the remainder of {q', I', t' , j\ A^^(A) |g, I, t, j)^^- 

After simplification, the complete contribution to {q' ,1' ,t' , j\ Ai'^{A) \q,l,t,j) from the instantaneous 
interaction below the cutoff is 



{q',l',t',j\MHA)\qJ,t,j)l 



— % / dx / dx' F(x,x'), 

167r4 A x-x' ' 



where 



F{x,x') = J 



= I d^k,d^k\—^Wn.„,e-^'"^"-'e-^'''"^^''^"'{x + x'){l-x + l~x'). 



X — x' 



To extract the divergence, we integrate Eq. ( |l5l| ) by parts with respect to x' : 

-2 /-l-e 



(g',r,t',j|X^(A)|g,/,Z,j)fj, 



167r4 



dx 



log(a; — x')F{x, x') I 



+ \og{x-x')F{x,x')\ + r ' dx'\og{x-x')^^^^p^ 

\x'=e dx 

= Bi + B2 + B^. 
The first contribution to {q' ,1' ,t' , j\ M'^{A) \q,l,t,j)f^ is 



167r4 e J2, 
X (1 - a; + i - 



dx I d^k^d^k'^iWq,q,e-^ e-^"" {x + x') 



(151) 
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(154) 
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To simplify this, we change variables from x to y: 

a; = y(l - 3e) + 2e, (155) 



and from kj_ and fc^ to Q± and N± 



— 2 ' 



Then 



X (2-2?/[l-3e]-3e)| ^_ . (157) 



I a; —x~e 



As e — ^ 0, the only contribution that survives is 



1 = T^^^'^'^'^S^l dyy{l~y)L\?\y)Ll^\y) J d'Q^4'\Q/2)Ti'\Q/2) J d^N^e-^^-'^' 

5,,,"^t,t'^A/|A2 C dxL\f\x)Lf{x)x{l ~ x)\oge. (158) 



Since (ffig2l -^^(A) |ffi52)sE is the part of (5^2! -^^(A) |5i52)se with the loge [see Eq. (|ll2|)], from 



Eq. (|129|) we see that {q' ,V ,t' ,j\M^{K)\q,l,t,j)'^^is]nst {q' ,1' ,t' ,j\M\K)\q,l,t,3)l^withthe [\ogx- 
11/12] factor replaced with — loge. This means that 

B, = -{q',l',t',j\M\K) \qj,t,j)^^. (159) 

Thus the divergence in {q' ,V ,t' , j\M'^{A) \q,l,t,j)fj^ cancels {q' ,1' ,t' , j\M'^{A) \q,l,t,j)^j^. 
The second contribution to {q' ,1' ,t' , j\ M'^{A) \q,l,t,j)^^ is 



B2 = ^/ dx\og(x-x') dkdk'd-/kk'e(k)e(k')0h)e(2TT~-f)^^W„„,e-^ 

StT-^ J2e J X-X' ' 

X e-2^"'^^^^^^(a; + x')(l -a; + 1 -a;')| • (160) 



-4 A 2 



To evaluate this, we change variables from fc' to s = fc'/v^- Then the leading term as e ^ is 

B2 = -e'^^-^^^cTt'^oAJo^y dx{2-x)logx J dkdsd-fkse{k)e{s)eij)e{2Tr-j)L'i^\x) 
X T/'''(/fc)W^5,,,e"''"[''+^'(^] 

= 0, (161) 
since e > —1/2. 

To simplify the third contribution to {q' ,1' ,t' , j\ A4'^{A) \q,l,t, j)fj^, we take the derivative of F{x,x') 
and take e ^ 0. Then 
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Nog: 

B3 = 
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f D\ogix-x')W,,,,T,^f\k')L['\x)T,^'\k)e-^-'^^-^^^ fl E„,, (162) 
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where 
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{El — dEi/dx', except for i = 2) and 



dL\?\x') 
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m'=0 

I' 



and 



dx' (1 - x'Y {x - x'Y 
dAiK _ (A?_i - fc|)2 

da;' a;'^ (2; — x'y 



This means that we can write the contribution to {q\l' ,t' ,j\M.'^{A)\q,l,t^j) from the instantaneous 
interaction below the cutoff as a divergent part and a finite part: 

{q' ,l\t' ,j\M\K)\q^,t,j)^^ = - {q' ,l\t\j\M\K)\q,l,t,j)'^^ + {q\l' ,t' ,j\M\K)\q^,t,3)^^ , (167) 
where 

{q'J'^',j\M\K)\q,l,t,j)'^i = j D\og{x-x')W,.,,,Tlf^{krLY\x)Ti''\k)e-''-''-^^^^^ 

5 5 

^ E^«' n (168) 

2—1 m— 1; mT^i 



Using the results of this section, the expression in Eq. (|126| ) for the matrix elements of the IMO 
becomes 



35 



(g', t',j\M\A) \q, I, t,j) = {q\ I', t',j\M\A) \q, I, t, j)ke + l'^ t'^j\ M\K) |g, 1, j)se 
+ {q', I', t',j\ M^{A) \q, I, t, j)cON + l'^ t'^ J\ M\A) \q, I, t, 

+ {q',l',t',j\M'{A)\q,l,t,j),^^^^ + {q',l',t',j\M^{A)\q,l,t,j}f^ . (169) 

Each of these terms is finite and we have taken e — > everywhere. We have written the first two terms as 
sums that can be computed numerically, and the four remaining terms as five-dimensional integrals that 
can be grouped into one integral suitable for numerical calculation. (See Appendix C for a discussion 
of some of the technical issues involved in the numerical calculation of these matrix elements.) Once 
we have computed these matrix elements, we can diagonalize the matrix to obtain glueball states and 
masses. 

6 Results and Error Analysis 

In this section we diagonalize the IMO matrix, obtaining glueball states and masses. We then discuss 
the sources of error in the calculation. We begin by discussing the procedure that we use to calculate the 
results. 

6.1 The Procedure for Calculating Results 

We represent a state using the notation Jf^ , where J is our best guess for the spin of the state, P is 
our best guess for the parity of the state, C represents the charge-conjugation eigenvalue of the state (it 
is always + because we have an even number of gluons), and j is the projection of the state's spin onto 
the 3-axis. We need to distinguish states with identical J's and P's and different j's because we do not 
have manifest rotational symmetry. If J = 0, we omit the subscript j in the state notation. We use an 
asterisk in the state notation next to the value of C to denote an excited state with the given quantum 
numbers. We will base our guesses for J and P on the numerical degeneracies of the states that have 
identical ,/'s and P's and different j's, and on the ordering of the states according to lattice data (see 
the discussion below). We consider only the five lightest glueballs (not counting as distinct those states 
that differ only in their value of j). The five lightest glueballs have spins J < 2. This means that we 
need to consider only \j\ < 2. For a given the states with j = \j\ and j = are degenerate and 
simply related (see Appendix C); so we explicitly consider only j = 0, 1,2. We consider nine values of 
the coupling: = g^/{4:iT) = 0.1, 0.2, 0.3, . . . , 0.9. To calculate our results, we implement the following 
four-step procedure. 

We execute the first step for all pairs (j, q,^). In this step, we define A = 1 and diagonalize 
{q',l',t',j\M^{A)\q,l,t,j) with all 4 spin basis functions (g = 1,2,3,4), but with only the lowest 
transverse-magnitude basis function {t = 0) and the two lowest longitudinal basis functions (/ = 0, 1). (It 
is necessary to use an even number of longitudinal basis functions so that symmetry and antisymmetry of 
the wave function under x ^ 1 — x are equally represented.) We perform this diagonalization as a func- 
tion of the basis-function parameters d and e that determine the widths of the transverse- magnitude and 
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longitudinal wave functions, respectively, and we find the values of d and e that minimize the ground-state 
mass. This yields what we consider to be the optimal wave- function widths for each pair (j, a^). 

We also execute the second step of the procedure for all pairs (j, a^). In this step, we fix d and e 
to be their computed optimal values and again define A = 1. We diagonalize the matrix with all 4 spin 
basis functions, A^t transverse-magnitude basis functions, and A^i = 2Nt longitudinal basis functions, for 
a total of 87V^^ basis functions, with A't = 1,2,3,..., 10. We use twice as many longitudinal functions 
as transverse-magnitude functions because \q,l,t,j) is zero if ^ + j is even and q = 4, or if Z + j is odd 
and q ^ 4. We want to use as many basis functions as possible, but we find for all pairs (j, a^) that 
when Nt > 7, the statistical errors from the Monte Carlo integrations of the matrix elements become 
overwhelming and the spectrum and wave functions become unreliable. The evidence of the breakdown 
is sudden contamination of the low-lying wave functions with high-order components. (See Appendix C 
for a more complete discussion of this topic.) Thus in the remainder of our procedure, we analyze the 
results that we find in this step when A't = 7, which corresponds to 392 basis functions. 

In the third step of the procedure, we use the mass of the ^ state (our most numerically reliable 
state) to determine the value of A for each . To do this, we note that A is the only mass scale in the 
problem. This means that the mass of the ^ state, Mq-+, can be written 

Mo-+ = 6(aJA, (170) 

where 6 is a dimensionless function of . Since we defined A = 1 in the second step of our procedure, 
the diagonahzation of the IMO as a function of yielded 6(a/^). In this third step, we consider A to be 
a parameter and define Mq-+ to be a constant. Then for a given coupling, we can use the results of the 
second step of our procedure to write the cutoff in units of 

A/A/0-+ = jj^. (171) 

Figure 5 shows the result for the third step of our procedure: a plot of the coupling as a function 
of the cutoff. When > 0.7 it is not a single- valued function of the cutoff. This is an indication that 
the coupling is too large. For this reason, we consider only < 0.7 in the remainder of our procedure. 
When < 0.7, the coupling decreases as the cutoff increases, as expected. However, depends on 
A more strongly than one may expect. We expect that perturbative pure-glue QCD would indicate 
that ~ 1/lnA, but the result that we get is much closer to ~ exp(— aA), where a is a constant. 
The reason for this is that the truncation of the perturbative series for Al^(A) and the truncation of 
the free-sector expansion of the states introduce spurious cutoff dependence in our results for physical 
quantities. Mq-+ is one such physical quantity. The spurious cutoff dependence of Mq~^+ is manifested 
through incorrect dependence of 6(q!a) on a^. This means that A has to compensate by depending on 

incorrectly in order to keep AfQ-+ a constant function of a^. This results in the strong dependence 
of on A that is shown in Figure 5. 

Using the recent anisotropic Euclidean lattice results of Morningstar and Peardon |^ , we can make 
a rough estimate of the range over which our cutoff is varying in Figure 5. They found that the mass of 
the ^ state is Mq-+ = 2.590 ± 0.040 ± 0.130 GeV. This means that our cutoff is varying from about 
3.1 - 6.0 GeV in Figure 5. 
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The fourth step of our procedure is to determine the optimal value of the cutoff, or equivalently, the 
optimal value of the coupling. We use two criteria to determine this. First we determine the value of the 
cutoff for which the computed masses are most independent of the cutoff. Figures 6, 7, and 8 show the 
masses of our states with j = 0, 1, 2, respectively, as functions of the cutoff. The masses and the cutoff 
are displayed in units of Mq-+. (Recall that Mq-+ was defined to be independent of the cutoff in the 
process of defining the cutoff as a function of the coupling.) The seven values of the cutoff at the points 
that we display correspond, from right to left, to =0.1, 0.2, 0.3, . . . , 0.7. ft is difficult to tell from these 
plots where the cutoff dependence is weakest (more points are needed) , but we see that the dependence 
is relatively weak from = 0.5 to ~ 0.7, which corresponds to A/Mq-+ = 1.33 to A/AIq-+ = 1.20. 

The second criterion that we use to determine the optimal cutoff is the degree to which the states with 
a given J and P and different j's are degenerate. This determines the cutoff that minimizes the violation 
of rotational symmetry. We find that these degeneracies are best when = 0.5. Given this and the fact 
that Figures 6-8 indicate that the cutoff dependence of the masses is weak when — 0.5 to 0.7, 
we determine that — 0.5 is the optimal coupling, and thus the optimal cutoff is A/Mq-+ = 1.33. 
Using the result of Morningstar and Peardon for the mass of the ^ state, we estimate that this cutoff 
is about 3.4 GeV. 



Now we present our main results. Our glueball masses for = 0.5 are summarized in Table 1, in units 
of the mass of the ground state (the 0++ state). Table 1 also shows an average of lattice results from 
a number of different calculations for the sake of comparison |^ . The uncertainties in our results that 
we report in Table 1 are only the statistical uncertainties associated with the Monte Carlo evaluation of 
the matrix elements of A^^(A). The full errors are much larger (see the discussion of sources of error 
below). We list three values of the masses for the 2++ and 2++* states for our calculation, corresponding 
to j — 0, 1,2. In each case the three masses would be degenerate if our calculation were exact. Our 
results agree with the lattice results quite well, perhaps better than we should have expected. 

We display our spectrum graphically in Figure 9. The masses are plotted in units of the mass of 
the 0++ state and the vertical widths of the levels represent the statistical uncertainties in the masses. 
The black lines connect the states that we believe should be degenerate. We see that the 2^^ and 2^"*"* 
glueballs are relatively degenerate with their 2^^ and 2^^* counterparts, and the 2^"''* is not too bad, 
but the 2^+ glueball is much too light. Our labeling of the states and subsequent assignment of the 
expected degeneracies are based on the ordering of the lattice states (see Table 1) and the apparent 
degeneracies of the 2^^ and 2q^* states with the 2^^ and 2^^* states, respectively. 

We want to show some of the features of the glueball wave functions. Rather than presenting the 
spin-dependent wave functions themselves, we present more illuminating spin-independent probability 
densities. A glueball state has the plane-wave normalization shown in Eq. (|8|) as long as the wave 
function ^i'^g^{x,k±) satisfies 



6.2 Results 




(172) 



This implies that 
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(173) 



where we define the dimensionless probability density 11 (a;, fc/A) by 

n{x,k/A)=27rAkY,\H:sA^,k^)\^- (174) 

SlS2 

We show the probability densities for some of our glueballs in Figures 10 - 14. The masses of the 
states tend to increase as the probability densities move away from the region a; 1/2 and more towards 
the edges. This is what we expect based on the form of the kinetic energy of a free state. Notice that 
the probability density for the 2 J"*" glueball is peaked around the region .t ^ 1/2 and looks similar to the 
probability density for the 0+"'" glueball. This is consistent with its small mass. 

6.3 Error Analysis 

We now turn to a discussion of the sources of error in our calculation. The sources of error are: 

• truncation of the renormalized IMO at 0{g^); 

• truncation of the frcc-scctor expansion of physical states at two gluons; 

• truncation of the basis-function expansion of wave functions; 

• numerical approximation of the matrix elements of Al^(A). 

Wc do not know how to estimate the size of any physical effects that require nonpcrturbative renor- 
malization. However, wc can naively estimate the size of the effects of higher-order pcrturbativc rcnor- 
malization. We have calculated the matrix elements of A^'^(A) through ©(o;^); so the corrections to 
these matrix elements should be 0(q;^). This translates to corrections to the mass spectrum of 0{a'^/2). 
Since we have used = 0.5, we estimate that the uncertainty in the mass spectrum from the effects of 
higher-order perturbative renormalization is about 13%. 

Wc do not know how to estimate the size of any physical effects that require an infinite number of 
particles. In fact, until wc include at least two free sectors, it is impossible to (iircc;tly estimate the size 
of corrections from higher free sectors. However, we can use the lack of degeneracy of the 2^^ state 
with the 2q and 2^"'' states to estimate these corrections. According to Table 1, the discrepancy in the 
various 2++ states is about 33%, if we believe the quoted lattice result. Since the uncertainty in the mass 
spectrum from effects of higher-order perturbative renormalization is around 13%, an uncertainty of 30% 
due to the truncation of the free-sector expansion is necessary to explain the lack of rotational symmetry 
in the spectrum (neglecting the other sources of error, which we expect to be small). 

As we mentioned, when we increase the number of basis functions that we use to represent the wave 
functions, we find that we reach a point where the statistical errors from the Monte Carlo evaluations 
of the matrix elements are overwhelming and cause our results to become completely unreliable (see 
Appendix C). For this reason, we have to truncate our basis-function expansion for the wave functions at 

= 7 transverse-magnitude functions. {N\ = 2Ni = 14, and there are 4 spin basis functions, for a total 
of 392 basis functions.) This truncation results in additional errors in our results. In Figures 15, 16, and 
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17, we show the convergence of the masses of the states with j — 0, 1, 2 respectively, as functions of iVt. 
The masses do not decrease as rapidly as functions of the number of states as one might expect. This 
is primarily because we have already optimized the states quite a bit by determining the widths of the 
transverse-magnitude and longitudinal basis functions that minimize the mass of the 0+"*" state (using 
A't = 1). Our best guess for the uncertainty that we introduce into the spectrum when we truncate the 
basis-function expansion, based on Figures 15 - 17, is a few percentf^. 

We can estimate the uncertainty in our results associated with the Monte Carlo evaluation of our 
matrix elements. To do this, we compute our results with a^^ =^ 0.5 four times, obtaining statistically 
independent results, and we compute the standard deviations of the masses that we obtain. This leads 
us to estimate that the uncertainty in the spectrum from the Monte Carlo routine is 1 - 2%. This is the 
uncertainty that we report in Table 1 and Figure 9. Combining this uncertainty with the others leads us 
to estimate that the total uncertainty in our results is about 33%. 



'^ Technically, this is not an uncertainty because improving the states can only reduce their masses, according to the 
variational principle. However, our discussion of errors is not meant to be rigorous. 
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Figure 5: The coupling as a function of the cutoff. We show the cutoff in units of the mass of the 
state, and we use 14 longitudinal basis functions, 7 transverse-magnitude basis functions, and 4 spin 
basis functions, for a total of 392 basis functions. Using the recent anisotropic Euclidean lattice result 
of Morningstar and Peardon for the mass of the state , we estimate that the cutoff is roughly 
varying from about 3.1 - 6.0 GeV in this figure. 
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Figure 6: The masses of the five lightest glueballs with j = 0, as functions of the cutoff. The masses 
and the cutoff arc displayed in units of the mass of the ^ state. The seven values of the cutoff at the 
points that wc display correspond, from right to left, to = 0.1, 0.2, 0.3, . . . , 0.7. We use 14 longitudinal 
basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, for a total of 392 basis 
functions. 
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Figure 7: The masses of the two lightest glueballs with j = 1, as functions of the cutoff. The masses 
and the cutoff arc displayed in units of the mass of the ^ state. The seven values of the cutoff at the 
points that wc display correspond, from right to left, to = 0.1, 0.2, 0.3, . . . , 0.7. We use 14 longitudinal 
basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, for a total of 392 basis 
functions. 
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Figure 8: The masses of the two lightest glueballs with j = 2, as functions of the cutoff. The masses 
and the cutoff arc displayed in units of the mass of the ^ state. The seven values of the cutoff at the 
points that wc display correspond, from right to left, to = 0.1, 0.2, 0.3, . . . , 0.7. We use 14 longitudinal 
basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, for a total of 392 basis 
functions. 
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State 


M/Mo+^ 


Lattice [28| 


0-+ 


1.38 ±0.02 


1.34±0.18 


2++ 


1.58 ±0.01 
1.58 ±0.02 
1.11 ±0.01 


1.42 ±0.06 


2++* 


1.70 ±0.01 
1.68 ±0.02 
1.62 ±0.02 


1.85 ±0.20 


0++* 


1.77 ±0.02 


1.78 ±0.12 



Table 1: The glueball masses from our calculation compared to an average of lattice results from a 
number of different calculations |2^. We display the masses in units of the mass of the O"*""*" state. 
The uncertainties for our results are only the statistical uncertainties associated with the Monte Carlo 
evaluation of the matrix elements of A4^{A). The three values of the masses for the 2++ and 2++* states 
for our calculation correspond to j = 0, 1, 2. We use the optimal coupling = 0.5, with 14 longitudinal 
basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, for a total of 392 basis 
functions. 
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Figure 9: Our glueball spectrum. The masses are plotted in units of the mass of the state and 
the vertical widths of the levels represent the statistical uncertainties in the masses. The black lines 
connect the states that we believe should be degenerate. We use the optimal coupling = 0.5, with 
14 longitudinal basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, for a 
total of 392 basis functions. 
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Figure 10: The probability density of the 0++ glueball. We use the optimal coupling = 0.5, with 
14 longitudinal basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, for a 
total of 392 basis functions. 
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Figure 11: The probability density of the glueball. We use the optimal coupHng = 0.5, with 
14 longitudinal basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, for a 
total of 392 basis functions. 
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Figure 12: The probability density of the 2o glueball. We use the optimal coupling = 0.5, with 
14 longitudinal basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, for a 
total of 392 basis functions. 
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Figure 13: The probability density of the 2^ glueball. We use the optimal coupling = 0.5, with 
14 longitudinal basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, for a 
total of 392 basis functions. 
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u.o ^ 

Figure 14: The probability density of the 22"*" glueball. We use the optimal coupling = 0.5, with 
14 longitudinal basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, for a 
total of 392 basis functions. 
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Figure 15: The masses of the five lightest glueballs with j = 0, in units of the cutoff, as functions of the 
number of transverse-magnitude basis fimctions, Nt. We use the optimal coupling = 0.5, with 4 spin 
basis functions and N\ = 2N^ longitudinal basis functions, for a total of SN^ basis functions. 
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Figure 16: The masses of the two hghtest glueballs with j = 1, in units of the cutoff, as functions of the 
number of transverse-magnitude basis hmctions, Nt. We use the optimal coupling = 0.5, with 4 spin 
basis functions and N\ = 2N^ longitudinal basis functions, for a total of SN^ basis functions. 
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Figure 17: The masses of the two hghtest glueballs with j = 2, in units of the cutoff, as functions of the 
number of transverse-magnitude basis hmctions, Nt. We use the optimal couphng = 0.5, with 4 spin 
basis functions and N\ = 2N^ longitudinal basis functions, for a total of SN^ basis functions. 
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7 Conclusion 



We have presented a formalism for pure-gluc QCD that aUows the physical states of the theory to rapidly 
converge in a free-sector expansion. In this approach, we force the free-state matrix elements of the IMO 
to satisfy three conditions to make the desired expansion possible. First, the diagonal matrix elements 
of the IMO must be dominated by the free part of the IMO. Second, the off-diagonal matrix elements of 
the IMO must quickly decrease as the difference of the free masses of the states increases. Third, the free 
mass of a free state must quickly increase as the number of particles in the state increases. 

We assume that we can use perturbation theory to derive the operators of the theory, and if this 
is valid, then the first condition is automatically satisfied. To satisfy the second condition, we place a 
smooth cutoff on the IMO. We use LFFT so that the effects of the vacuum are isolated in particles with 
zero longitudinal momentum, and we remove these particles from the theory with the intent of replacing 
their physical effects with interactions. This makes it reasonable to expect that the third condition on 
the IMO is satisfied automatically due to the free-particle dispersion relation of LFFT. 

The cutoff that we use violates a number of physical principles of light-front pure-glue QCD. However, 
by requiring the IMO to produce cutoff-independent physical quantities and by requiring it to respect 
the unviolated physical principles of the theory, we are able to derive recursion relations that uniquely 
determine the IMO to all orders in perturbation theory. 

We have applied our method to the calculation of physical states and masses. For this calculation, 
we approximated all physical states as two-gluon states. We calculated the color parts of the states 
analytically, and we expanded the states' momentum and spin degrees of freedom in terms of basis 
functions. We designed the states to be simultaneous eigenstates of the IMO, the three-momentum 
operator, and the projection of the internal rotation generator onto the 3-axis. 

Using our recursion relations for the IMO, we calculated to second order in perturbation theory the 
two-gluon to two-gluon matrix element of the IMO, which is required for the calculation of physical states. 
We then used it to calculate the IMO matrix in terms of the basis functions. We showed that the infrared 
divergences in the matrix from exchanged gluons with infinitesimal longitudinal momentum cancel when 
treated properly. 

In order to diagonalize the IMO matrix, we computed the five- dimensional integrals in the matrix 
elements rising Monte Carlo methods. We calculated the ghieball spectrum for a range of couplings and 
fomid that we could not use more than about 400 basis functions without the statistical errors becoming 
overwhelming. We used the mass of our ^ glueball to compute the nonperturbative cutoff dependence 
of the coupling, and we analyzed the cutoff dependence of the spectrum. We found that the cutoff 
that minimizes our errors is A/Mo-+ = 1.33. The corresponding coupling is a,^ = 0.5. We presented 
the probability densities for some of our glueballs and found that our results for the spectrum compare 
favorably with recent lattice data. The largest discrepancy seems to be the state, which is much too 
light. Finally, we analyzed the errors in our calculation from the various possible sources, and estimated 
the total imcertainty in our spectrum to be 33%. 

There are two main paths that we can take for future work with our approach. The first path is 
to further test our method with the theories that we have considered so far. Since the scalar theory 
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that we considered in Ref. |l| is relatively simple, it would be interesting to use it compute the IMO 
to higher orders in perturbation theory. This would require us to solve the integral equation for the 
cutoff-independent reduced interaction and could be used to further check our conjecture that our IMO 
leads to correct scattering amplitudes order-by-order in perturbation theory. 

In pure-glue QCD, we can further test our approach for calculating physical states by computing the 
IMO to higher orders in perturbation theory and by keeping more free sectors in the expansion of the 
states. However, to keep more free sectors in the expansion of the states, we have to be able to calculate 
IMO matrices that have more degrees of freedom, while controlling the statistical errors in the spectrum. 
This means that we need a better algorithm for determining how accurately individual matrix elements 
of the IMO have to be computed in order to get a desired uncertainty in the spectrum. We could also use 
a better basis that requires fewer momentum functions to represent a wave function. Overcoming these 
problems will be challenging, but it is important to test our method by studying the rate of convergence 
of the free-sector expansion of states as a function of both the cutoff and the masses of the states. 

Another test of our method that we can do with pure-glue QCD is to analyze the interaction that 
we have derived in this paper to test the conjecture that it is logarithmically confining. It would also be 
interesting to analyze the long-range parts of higher-order interactions to see if the perturbative series 
for the interaction is building towards a linearly confining potential. Analyzing the long-range parts of 
higher-order interactions may be much easier than computing these interactions in their entirety. 

The second main path that we can take in the future is to extend our method to other theories and 
operators. In order to compute quantities that can be compared with experiment, we wish to extend 



and numerical complexity from the the vertices involving quarks. Second, quark masses complicate the 
method for determining the IMO because they increase the number of reduced interactions that can be 
cutoff-independent jist . In addition, if large and small quark masses are considered simultaneously, then 
efficient numerical representation of the states and accurate calculation of the IMO's matrix elements 
become more difficult. Masses also quickly enlarge the parameter space that must be explored to compare 
to experimental data. 

We can also extend our method by applying it to the computation of operators other than the IMO, 
such as the rotation generators, the parity operator, and currents. The rotation generators and the parity 
operator are of particular interest because they may aid in the classification of the physical states of a 
theory. 

In summary, there are many avenues of research that must be explored, and some of them are quite 
complex. However, all the improvements that we have discussed are necessary if we are to accurately 
represent the physical states of quantum field theories as rapidly convergent expansions in free sectors. 



^We are not thinking of including QCD eflfects that require nonperturbative renormaUzation or an infinite number of 
particles because these would require a new method rather than an extension of our current approach. 



our method 




complicated for two reasons. First, there is additional algebraic 
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APPENDIX A: Conventions for Light-Front Pure-Glue QCD 

The purpose of this appendix is to state our conventions. With these conventions, any four-vector a 
is written in the form 

a = (a^, a^, aj_), (175) 

where in terms of equal-time vector components, 

a^=aP±a?, (176) 

and 

2 2 



where is the unit vector pointing along the z-axis. The inner product is 

a • = + ia"5+ - a^- (178) 

and 

d\ = a^^ - d^. (179) 
A spacetime coordinate is a four- vector, and according to Eq. ( |l75| ), it is written 

x = {x+,x-,xi_). (180) 

The time component is chosen to be x'^ . x~ is referred to as the longitudinal component, and x±_ contains 
the transverse components. 

The gradient operator is treated just like any other four-vector. Its components are 

^^ = 2^> (181) 



and 



The canonical Lagrangian density for pure-glue QCD is 

^ - -^^c^.^r, (183) 

where 

= a^^c - c'^A^ - 5^^, (184) 

Greek indices are Lorentz indices, c's are color indices, repeated indices are summed over, and the /'s 
are the SU(A'c) structure constants. 

We derive the canonical Hamiltonian from C by the following procedure: 

1. We choose the light-cone gauge, — 0. 
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2. We derive the Euler-Lagrange equation that determines in terms of A±c. 

3. Using the canonical procedure and treating the field classically (i.e. letting it and its derivatives 
commute), we derive the Hamiltonian in terms of A±c, dropping terms that are zero if the gluon 
field is zero at spacetime infinity. 

4. We quantize the gluon field by expanding it in terms of free-particle creation and annihilation 
operators. (We define the field expansion and its inverse longitudinal derivative below.) 

5. In each term in the Hamiltonian, we treat the creation and annihilation operators as if they commute 
and move all the creation operators to the left of all the annihilation operators. (This "normal 
ordering" drops the so-called "self-inertias," as well as some constants.) 

6. We drop the terms in the Hamiltonian that have no effect if there are no particles with j3+ = 0. 

Wc work in the Schrodinger representation, where operators are time-independent and states are 
timc-dcpcndcnt. Thus we quantize the field by defining it to be a superposition of solutions to the 
Klein-Gordon equation (since gluons are bosons), with the quantization surface a;"*" = 0: 



(185) 

x+=0 



where 

Ci=lsi=-l,l 

Here Si is the spin polarization of particle i, P is the four-momentum operator, e is a positive infinitesimal, 
and the gluon polarization vector is defined by 

£1, = (187) 

p+ and p_L are the momenta conjugate to x~ and x±; so they are referred to as the longitudinal and 
transverse momenta, respectively. The purpose of e is to regulate divergent effects from exchanged 
gluons (either instantaneous or real) with infinitesimal longitudinal momentum. We take e — > before 

we calculate physical quantities (sec Sciction 5). 

In the process of deriving the canonical Hamiltonian, wc need to take the inverse longitudinal derivative 
of the gluon field. We do this by moving the derivative inside the expansion: 



-^A±c{x ,x±) = I DiSc 



(188) 

x+=0 



The expansion coefficients a| and Oj are identified as particle creation and annihilation operators. 
They follow the convention 

ai = a{pi,Si,Ci), (189) 

and have the commutation relations 

[ai,a]] = Si J = 16Tr^p+ 5^^^ {pi-pj)6s,,sjSc„cj (190) 
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and 

[a„aj]^[ala]]=0, (191) 

where 

5(3) (p. -p.) = s{pt-p+)S^^\pi^-pj^). (192) 

Let A4^ be the invariant-mass operator. The momentum conjugate to a;"*" is p~; so the Hamiltonian 
is identified as 'P" and it follows from 

= M"^ (193) 

that 

V- = ^^^r^- (194) 

The canonical Hamiltonian that results from the above procedure can be written as the sum of a free 
part and an interacting part: 

where 

^f7ee = / D,?^aW (196) 

J Pi 

We give the interaction in terms of its "modified" matrix elements: 

V = J D1D2D3 44ai 5(3) (pi -p2- Pa) ^;S'25'3|^^|5'i^ 

+ 9^^- J D1D2D3 alaia2 S'-^'^ {pi +P2- Ps) ^93^9192:^- 

+ 9^'§^ I DiD2D3D4aWiaia2 S^^\pi + P2 - Ps - Pi) ^g3g4|^^|gig2:|>» 

' ' ^ i=l 

+ f^-^- / DxD2D3DAalala\ai5^^\pi-p2-p3-PA)^^92g3gA\v\gi'^i 

■ 1=1 

+ 9^^^ j DiD2DsDia\aia2as5^^\px+p2+P3-pA)^^9M9i9293'ku (197) 

■ i=l 

where a modified matrix element is defined by 

(198) 



167r35(3)(ft - ft ) 
The modified matrix elements are 

^9293\v\9i> = ir^'^' 
+ S. 



3=1 



1,83^182- |(Pl-L +P3-L) - ^(Pf +^^3 )| 



P2 
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I P3 



Ss2,S3^±si ■ \ {P3± +P2±)- ^{Pt +P2) 

I Pi 

+ 5s,,s3S±S2 - ^-{Pl± + P3±) + ^{Pl + Pt)^ 



(199) 



and 



ta4i^i5i52^2 = /-'=-r--^«„.3^..,./^^^^^''"^''''' 

^5354^15152:^3 = r'''T'"'Ss„sJs 



Pt 


-pt pt 


-Pt 


Pi +pt pt +pt 


Pt 


-ptpt 


-Pt 


Pt 


-pt pt 


-Pt 



^5354|^;|5i52^4 = r^-^r^^^^(?....-.^.3..-4 + . 1- . (200) 



'''^Pi +Pt Pi +Pt' 



and 



4:525354|t^i5i^i = ir^'^^'^r^"-" {-s,,,,,s,„s, + s,,,sjs,,s2) + r^'^^^r^"-" {-s,,,sjsus. + Ss„sjs„s2) 

, fClC2CfC3C4CfSj Sj , Sj S- )] 

^525354^151^2 = -r^^^7^^^^^'5,„,,(5,3,,,4±44^' 

Pi - P2 Pi + Pt 

^525354|t^|5i^3 = -r^^^V^^^^^^.„«3^«-.,.4 4^4^' 

Pi -P3 Pi + P2 

i525354|t^|5i^4 = -r'''r''''Ss,,sJ,3,s2 ^i^^i , (201) 

Pt-PiP2+Pt 

and 

^541^^1515253^1 = [.n'''n''^-6s2,,36s,,a,+6s,,sJs,,S2) + n"'" r'"'" {-5s2,sJs,,S. + 5s3,sJs.,-S2) 
I fClC2C fC3C4C X _ S - )] 

^54|t^|5l5253:^2 = -r^^*V'^'^'^ai,«4^a 

4:54^1515253^3 = -r''''r''''S,,,,J, 

^54|t^|5i5253^4 = -/c3C4cyc.c2c^^^^^^^^^ _^ P3 ^''X'\~/l , (202) 

P3 -PiPl +P2 

where s„ = — s„. 

We work in the free basis, the basis of eigenstates of V^^^. They are given by 

|5i52---5n) =44---4|0), (203) 



pt + pt pt 


-Pt 


pt - pt pt 


+ pt 


pt + pt pt 


-Pt 


pt - Pi pt 


+ pt 


pt + Pi Pi 


-Pt 
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for any integer n > 0. The associated eigenvalue equation is 

^froc 15152 ■■■9n) ^9192 ■ ■ ■ 9n) , (204) 

where 



PiA 



p- = ^ (205) 



(since pf — 0) and the sum is zero if n = 0. 



The noninteracting hmit of Eq. (194) is 

Pf+_Ml 
V 



n;ee- ^ (206) 



where Mf^^^, is the free invariant-mass operator. It has the eigenvalue equation 

Ml,, |5i<?2 • • • g„) = 15132 • • • gn) , (207) 

where 

M^ = P+"Ypr -Pl, (208) 

and P is the total momentum of the state. 

Finally, in terms of the free states, the completeness relation is 

1 = |0) (0| + I Di (5i| + ^ / D,D2 I5152) (3152! + • • • . (209) 
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APPENDIX B: The Derivation of the Recursion Relations for the Reduced 
Interaction 



This appendix is an extension of Appendix D of Ref. to the case of pure-glue QCD. In Subsection 
2.4 , we derived a constraint on the 0{g^) reduced interaction for r > 1: 



(210) 



Since we aheady know the first-order reduced interaction [see Eq. (18)] and the cutoff-independent part 
of the second-order reduced interaction [see Eq. (p^)], we wish to use this equation to compute Vqq(A) 

(r) . 

for r > 2 and VJAj for r > 3, in terms of lower-order reduced interactions. 

B.l The Recursion Relation for the Cutoff-Dependent Part 

We begin by computing the recursion relation for the cutoff-dependent part of the reduced interaction. 
Recall that momentum conservation implies that any matrix element of V^(A) can be written as an 
expansion in unique products of momentum-conserving delta functions. This means that an arbitrary 



matrix element of Eq. (210) can be expanded in products of delta functions: 



s=2 



-(211) 



where the (z) superscripts denote that we are considering the i^^ product of delta functions that can 
occur in a delta-function expansion of {F\ ^(''^(A) |J). This equation is equivalent to a set of equations, 
one for each possible product of delta functions: 



(J^|yW(A)|^)W_(^|yW(A')|j)W 



(F| <5T/('-) [/) - Y Br-s,s {F\ V^--'\A) \I) 



(i) 



Cluster decomposition implies that we can write 



(F|yW(A)|/)« 



F^%{{pn}ASn}ACn},A), 



(212) 



(213) 



where Sj is the j momentum-conserving delta function in the i product of delta functions (S 
includes a longitudinal- momentum factor), ivj*-* is the number of delta functions in the i^^ product, and 
Fq^^ is a function of the cutoff and the quantum numbers of the particles in the matrix element, but does 
not contain delta functions that fix momenta. We define A'part to be the number of particles in state [/) 
plus the number of particles in state \F), and n = 1, 2, 3, . . . , Apart- The momentum, spin polarization, 
and color of particle n are given by p„, s„, and c„. We define N^^^ to be the number of particles in the 
matrix element that participate in an interaction for the i*'' product of delta functions. In order for the 
IMO to have the dimensions (mass)^, must have the dimensions (mass)''~^in' . Note that we are 
suppressing the dependence of the RHS of this equation on r. 

We have assumed that any matrix element of the IMO can be expanded in powers of transverse 
momenta, not including the momentum-conserving delta functions; so 
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4D({Pn},{s«},{c„},A)=A4-< ^ ^ 

{m„t} 



{p+},K},K}) n n( Y 



Npart 2 / f \ rrint 
PnJ 



(214) 



where t denotes a component of transverse momentum and m„t is a non-negative integer index associated 
with transverse-momentum component t of particle n. The sum is over all values of each of the ninths, 
subject to the constraint that 



(215) 



which is necessary to avoid terms in the momentum expansion that are cutoff-independent. The 2;|™"''^'s 
are the coefficients for the momentum expansion. They depend on i and the ninths and are functions of 
the particles' longitudinal momenta, spin polarizations, and colors. 



Since the RHS of Eq. (212) has the same product of delta functions as the LHS, we can write 



(F| |/) - J2 Br-s^s {F\ y('-^)(A) |/) 









n 




i=i 



G«({p„},K},{c4,A,A'), (216) 



where G^*-' has dimensions (mass)'* ^int , and by inspection of the LHS of Eq. (|212D and Eq. (21S), is a 



function of the quantum numbers of the particles. A, and A'. Substitution of Eqs. ( ^16[) and (|213| ) into 
Eq. (|l|) yields 

F^{{Pn}, {sn}, {c„}, A) - F^liipr,}, {s„}, {c„}. A') = G« ({p4, {s„}, {c„}, A, A'), (217) 

where the momenta in this equation are constrained by the delta- function conditions. 

Since the LHS of Eq. (^17|) is the difference of a function of A and the same function with A A', 



G^*' must be as well. Since the LHS of Eq. (217) can be expanded in powers of transverse momenta, G*^*-' 
must have the form 



G«({Pn},{s„},{c„},A,A') = 4""'^({P^},{U,{cJ 

{m,it} 



A- 



n 



A 



A' 



n 



A' 



(218) 



where the sum is restricted by Eq. (215). 

Substituting Eqs. (|l|) and (|214| ) into Eq. (|l|), we find that 



A' 



{p+},K},K})nf47 



{m„t} 



{m„t} 



A' 



n 



A 



A' 



n 



Pnl 
A' 



(219) 
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Matching powers of transverse momenta on both sides of this equation gives 



A-^-^.-t-^"-' ™"* - A'4-<'t-s.,t m™tl , (220) 

The factor in brackets cannot be zero because A 7^ A' and Eq. ( |215| ) holds. Thus Eq. (220) impUes that 

(221) 



t"'\{pt}ASn},{Cn} 



Then Eqs. ( |214| ), ( pi8| ), and ( |22l| ) imply that 

4d({p4, {^^n}, {c„}, A) = G(') ({p„}, {s„}, {c„}, A, A') I 



A terms 



(222) 



where "A terms" means that G*-*-' is to be expanded in powers of transverse momenta and only the terms 
in the expansion that are proportional to powers or inverse powers of A contribute. From Eqs. ( |213| ), 
(IH), and (HI), 

r-l 1 W 

{F\ 5V^^^ ^^-^'^ (^1 ^^'■"'^ (^) 1^) ' (223) 

■^—2 J A terms 

where it is understood that the momentum-conserving delta functions are ignored for the purposes of 
transverse-momentum expansions. Since a matrix element is the sum of the contributions to it from 
different products of delta functions, both sides of this equation can be summed over i to obtain 



{F\V^^^{A)\I) 



{F\V^'^{A)\I) = 



{F\ |/) - J2 Br-s,s {F\ V^^-^\A) \I) 



(224) 



A terms 



This equation tells us how to calculate the cutoff-dependent part of the 0(5^) reduced interaction in 
terms of lower-order contributions. 



B.2 The Recursion Relation for the Cutoff-Independent Part 

Since we have specified V^-^^ and Vq^ , we need to determine Vq^ for r > 3. It is useful to first consider 
which contributions to V^'^''^(A) can be cutoff-independent. 

A matrix element of the cutoff-independent part of 1/'^'"^(A) can be expanded in products of delta 
functions and in powers of transverse momenta just as was done for the cutoff-dependent part. Thus we 
can write 

{F\V^\\l)^Y.^P\<^\^)''^^ (225) 

i 

where 

F^Gl{{Pn}ASn}ACn}) (226) 

and 

F^^l{{p^}.{s..}.{c.r^})=K'-^'^^^ Y: wl"^-'\{pt},{^n},{Cn}) J] 11 [^) ' (227) 

{m„t} n=l t=l ^ ^ 



{F\V^[^\I)^'^ 



n4 
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The sum is over all values of each of the nint's, subject to the constraint that 



(228) 



which ensures that all the terms in the expansion of pj^j are cutoff-independent. 



Eq. (228) places constraints on the possible cutoff-independent contributions to the reduced interac- 



tion. Any contribution to a matrix element of V"'''^(A) has an N^^^ > 2, but Eq. (228) can only hold if 



Suppose that N^^^ — 2- In this case, Eq. (228) implies that F^l is quadratic in transverse momenta. 
Due to approximate cluster decomposition, only interacting particles' transverse momenta can appear in 
FqI- So any contribution to F^j can depend on the transverse momenta of two interacting particles. 
Thus F(^{ can be written as a sum of terms, where each term corresponds to a distinct pair of interacting 
particles. The momentum dependence of each term in f^'j is limited to dependence on the momenta of 
the interacting particles and the total longitudinal momentum: 



^F^r^(fc„.,P+,{s„},{c„}), 



(229) 



where km and fc^ are the momenta for the initial and final particles in the m^^ interacting pair, and 
where we have used the fact that for N^^^ — 2, momentum conservation implies that km = k'^. F^^^ 
must be quadratic in fc,,,^ or it must be zero. The matrix elements of the IMO are boost-invariant, as is 
the delta- function product in Eq. ( ^26[ ). This means that F^l must be boost-invariant, but it cannot be 



if Fq^""' is quadratic in fcm±; so F^^""' must be zero. Thus the reduced interaction does not contain any 
cutoff-independent two-point interactions. 

Note that two-point interactions are self-energies, and they change the particle dispersion relation. 
If they change the dispersion relation such that the coefficients of the free relation become modified 
by interactions, then this can be viewed as renormalization of the field operators, i.e. wave-function 
renormalization. This effect is absent unless either F^\ or can be quadratic in transverse momenta 
for N^^^ — 2. We have just shown that boost invariance prevents this for F^l , and according to Eq. ( ^ 



cannot be quadratic in transverse momenta for A^V^ = 2; so there is no wave- function renormalization 
at any order in in our approach. 

According to Eq. ( |22^ ), if A^j^'^ — 3, then F^\ has to be linear in transverse momenta, and if A^j'^^^ = 4, 
then FqI has to be independent of all transverse momenta. According to assumptions that we made in 



(r) \ 

Subsection |2.3.4 if r is odd then Vq-^ has no A^j 



4 part, and if r is even then V^^j has no N^^^ 



(r) 

part. This means that if r is odd, then V^j can contain only three-point interactions that are linear 

(r) 

in transverse momenta, and if r is even, then V^^' can contain only four-point interactions that are 
independent of all transverse momenta. 

-r + 2: 



To calculate V^\\ we consider Eq. (212) with 



{F\V^^+'\A)\lf^ - {F\V^+'\a')\I) 



(0 _ 



r+l 



+2-S 



(A)|/) 



.(230) 
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(i) 

In the remainder of this appendix, we assume that r is odd. Then we need to consider only N-^^^ — 3 



initiaUy. We expand Eq. ( 230 ) in powers of transverse momenta and keep only the linear term: 



= 



r+l 



(i) 



(231) 



s—2 J term 

where we have used the fact that {F\ V^'j^'^\A) |/)^'^ has no linear part when N^^^ = 3 [see Eq. (pTs])]. 
We move the first term in the sum on the RHS to the LHS: 



BrAF\V^^\I)^ = 



r+l 



{F\ <5F(''+2) |/) - J2 Br+2-s,s {F\ ^(-+2--) (A) |/) 



s=3 



(i) 



pj- term. 



(232) 



Now we can sum over all values of i corresponding to three-point interactions: 

r+l 



{F\V^'^\I) 



1 



Br 



[F\ |/> - Br+2-s,s {F\ ^('■+2-^) (A) |/) 



3 — point 



(233) 



To use this equation, we also need to use Eq. ( ^30| ) to solve for V^j^^^. Since r is odd, r + 1 is even. 
Thus V'qj^^'' will contain only transverse-momentum-independent four-point interactions. Using steps 
analogous to those that led to Eq. (|233|), we find that 



Br 



-1-1,2 



r+2 1 4-point 

(F| 6V('-+'^ 1^) - E Br+3-s.s {F\ (A) |/) 

s=3 



(234) 



To use these equations, the right-hand sides have to be expanded in powers of transverse momenta. Only 

(r) 

three-point interactions that are linear in transverse momenta contribute to V^j , and only four-point 
interactions that are independent of all transverse momenta contribute to V^\'^^\ 

These equations are coupled integral equations^ because both Vq^^ and V^i^^^ appear on the RHS 
of Eq. (HI) inside integrals in 6V'^''+'^\ and V^\^^^ appears on the RHS of Eq. ( |234[) inside integrals in 
5y(''+3). It would seem that Vqj^^^ also appears on the RHS of Eq. (|234|) inside integrals in SV'^'''^^\ but 



^ci^^^ cannot couple to T^'^^ to produce a transverse-momentum-independent four-point contribution to 
j\/(''+3) This is because the cutoff function T2^'^ ■* vanishes when the intermediate state is massless and 
all external transverse momenta are zero. This means that since we specified t^*-^' and V^j^ in Subsection 
2.3.4 , we can use Eqs. (|233| ) and (234) to solve for V^j^ and V^j'* simultaneously, and Vqj^ and 
simultaneously, and so on. Note that before we can use these equations to solve for VJ^j-* and Vjlj^^-* 



simultaneously, we must first use Eq. (224) both to compute VqjI{A) in terms of lower-order interactions 
and to express Vq'q^^'(A) in terms of lower-order interactions and y('')(A). 

Before concluding this appendix, we would like to deduce a bit more about the relationship of 5^ to 
g^, . According to Eq. ( p7| ) and the surrounding discussion, this relationship is determined by the matrix 
element (52<73| 6V \gi), which can be expanded in powers of (7^,: 



(.9253 1 '^F \9i) =Y.9l (5253 1 '5^^^*^ \9i) 



(235) 



t=3 



'^^It is very difficult to prove that integral equations of this type have a unique solution; so we simply assume that it is 
true in this case. 
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Recall that SV'^*^ is built from products of t^^''^ (A')'s. This implies that SV'-*'^ can change particle number 
by 1 only if t is odd, and thus Eq. ( |235| ) implies that the coupling runs at odd orders; i.e. C's is zero if s 
is even [see Eq. (^]. 
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APPENDIX C: Technical Issues in the Numerical Calculation of Matrix 
Elements 



In this appendix, we discuss some of the technical issues involved in the numerical calculation of the 
matrix elements {q' ,1' ,t' , j\M.'^{K) \q,l,t,j). In Section C.l we discuss how we put the matrix elements 
into a form that is amenable to numerical calculation. In Section C.2 we briefly cover three topics: we 
show how the glueball state with J';^ eigenvalue —j can be written in terms of the glueball state with 
eigenvalue j, we list a few tricks that allow us to reduce the number of matrix elements that we must 
compute, and we present our method for estimating how numerical uncertainties in the matrix elements 
translate into uncertainties in the spectrum. 

C.l Preparation of Integrals for Monte Carlo 

There are two types of contributions to {q' ,1' ,t' ,j\M^{A) \q,l,t,j): finite sums and five-dimensional 
integrals. We use Mathematica to evaluate the finite sums to as many digits as we wish. To evaluate 
the integrals, we combine them into one integral and use the VEGAS Monte Carlo routine It takes 
a bit of work to put the integral into a form that will converge. 

There are two main difficulties with getting the integral to converge. The first difficulty is that 
{q' , I', t' ,j\M^{A) \q, I, t, looks divergent: as x' x, the sum in Eq. ( |16S| ) diverges. This divergence 
is misleading because it actually contributes nothing to the integral (assuming that we calculate the 
integral carefully; see the discussion below). Left unchecked, this false divergence prevents the integral 
from converging with VEGAS. To rectify the problem, we want to subtract the false divergence from the 
integrand. Since it integrates to zero, this is allowed. 

The second main difficulty with getting the integral to converge is roundoff error. Even after we 
subtract the false divergence from the integrand, the integrand peaks around x' — x, and in this region 
there are large cancellations in some of the quantities that we have defined. To prevent these cancellations 
from causing roundoff error, we rewrite these quantities so that the cancellations are explicit, before we 
turn to numerics. 

C.1.1 Subtraction of the False Divergence 

We begin by defining a set of variables that is natural for dealing with the false divergence. We define 

ri^x-x'. (236) 

We change variables from k±_ and to the dimcnsionless transverse variables rj_ and w± : 

r± = ^(^-i + fcl), 

We define the angle between r± and w± to be f3: 

r± ■ w± = rw cos /3. (238) 

We also define 
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r± = y +r]w^ ±2rwy/rj cos 13, (239) 
and then we can derive a host of useful relations: 



k 
k' 



d 
d 

!± 
d ' 

r_ 

T' 
1 



cos 7 



9 9 

r — rjw 



sin 7 = ^™^ sin/3, (240) 



and 



dkdk' djkk' e{k)0{k' )e{-f) 9 {2Tr -j) = ^ dr dw d(3r w 0{r) 0{w) 0{f3) 0{2Tr - P). (241) 

In terms of these variables, (q', t\j\ Al^(A) |g, i, j)f^ takes the form 

(g',r,t',j|A<2(A)|q,/,t,j)^^^ = y'i?'log^T4^,^,,r,,(r„)T,(r+)e-(^^)"'(^-+^-) 

5 5 

E^' n (242) 

2—1 m— 1; mT^z 

where 

D' = dx dx' drdwdPrwT] 0{x) 0{1 - x) 0{x') e{x - x') e{r) e{w) e{P) 6(271 - /3)Z|''^ (x). (243) 

To avoid roundoff error, we have defined simplified dimensionless versions of the differences of the free 
masses and the derivatives of these differences: 

r^r]{l — X — x') + w'^r]^{l ~ x ~ x') — 2wr^{x[l — x] + x'[\ — x']) cos (3 



Afi 
Afk 

AiK 



x{l — x)x' {\ — x') 
r'^rj + w'^{l — X + 1 — x'Y + 2rwy^{l ~ x + 1 — x') cos/3 
(l-x)(l-a;') ' 
r'^rj + w^{x + x'Y ^ 2rwy/rj{x + x') cos (5 



A'pj, - ^-4 



(1 — x')^ rj 
up 

^'iK = -71 -4-. (244) 



Eq. ( |242| ) is now dimensionless, except for the factor in front of the integral, which is proportional to 
l/d\ 
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As 



X (i] 0), the contributions to the integrand from the first and second terms in the sum 



diverge, fn the hmit x' -~> x, the contribution to the integral from these terms can be written^^ 



X (1 - a; + 1 - a;')- 



5,.,,, / i?'log7?T,,(r)T,(r)e-32(A'^)-^-^ 



4 - 64(Ad)-4!^ 



l'{;\x'){x + x') 



(245) 



An examination of this integral reveals a problem: the transverse integrals are zero and the longitudinal 
integrals are infinite. To solve this problem, we consider what would have happened if we had not yet 
taken e ^ 0. In this case, the transverse integrals would be zero and the longitudinal integrals would be 
finite. Thus the e ^ limit of this integral would be zero. This means that this integral is actually zero, 
and we can subtract it from the full integral in Eq. ( ^42| ): 



q',l\t',j\M\A)\q,l,t,j)f^ 



D' log T] 



'^k+^'ik) J2 El Yl E„ 

2—1 m— 1; m^i 



1 7/1^ 

4-64(Ad)-4!^ 



L[?\x')ix + x')il-x + l-x') 



(246) 



Once we have performed this subtraction, there is no ambiguity about the value of the full integral, and 
it converges when computed numerically. 

C.1.2 Combination of the Integrals 

We now combine all the five-dimensional integrals into one integral: 



{q',l\t\j\M'{A)\qJ,t,j),_^ 



27r3d2 



D' 



Ic 



ON + ^EX + ^IN+EX 



'IN 



(247) 



where 



CON 



L|f' (x')T.(-)T^*(r+)e-(^'^)-^^-i f + A^) (l 



_ g-2(A<i)-''Af K 



+ 



'IN+EX 



X{\-X) «•« X'(l-X') X{\- X)x'{\~ X') 

"A^ V 

^FI _1 



' ^ ^ ^ / n +7 x{\ - x)x' {\ ~ x') r]AFK^iK 

- rfw^ - 2Tj*w^{r^ + 2w'^[l - 2x]) - rf{r^ + Ar'^w^[l - 2x] + 6u;^[l - 2xf) 
4cos/3^rw([-l + rf]r'^ + r77«^[l - x - x'Y){\ - x - x') + 8r^w^x(l - x) 
Acos^lSr^w'^iri^ + 2r]^[l - 2x] +4a;[l - x] - 4r]'^x[l - x] ~ 2r][l - 2x]) 
iri^w^iw'^ll - 2xf +r^[l- 2x{l - x)]) 
2-q{r^ - 2r'^v?[l - 2x\ - Aw^x[-l + ix - Ax"^ + 2x^])\ , 



'We do not replace all the occurrences of x' with x because doing so hampers the convergence of the integral. 
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'IN 



i—1 m— 1; m^i 



log?7 



4 - 64(Ad)-4^ 



Z[f)(x')(x + x') 



X (1 - a; + 1 - .t') 



(248) 



We have rewritten the integrand of {q' ,1' ,t' , j\ A4'^{A) I?, ^ j)iN_|_Ex ehminate large roundoff errors, 
at the expense of making it more compUcated. To further avoid roundoff errors, we rewrite a few of the 



s. 



(1) 
3,3 

^(3) 



cos(j7)(2?72 + 2?7(1 - 2x) + (1 - 2xf), 



S^ i = -sin7sin(j7)(2r7^[l-2x] +77[l-6x + 6a;-^] -2a;[l-3a; + 2a;^]), 



'4,4 



cos(j7)(l - 2x)(l + 2t]- 2x). 



(249) 



Note that to compute some of the trigonometric functions that appear in this integral [such as cos(j7)] 
in terms of the integration variables, it is necessary to use recursion relations that define these functions 
in terms of cos 7 and sin 7 so that we can use Eq. (|240| ). 



The integral in Eq. (247) converges slowly. This is because it is strongly peaked when x' ~ x, even 
though we have subtracted the false divergence. To spread out this region, we change variables from x' 
to p where 



X = X 



(1 - e-) 



(250) 



Now the integral is 

{q',l\t\j\M'{A)\q,l,t,j),_^ = - 



2ti^<P 



dx 



dp 



dr 



dw I dp r w if L^f^x) 



IcO^ + IeX + -^IN+EX + ^i^^ 



(251) 



As a final step, we note that VEGAS requires the region of integration to be finite. Thus we change 
variables from p, r, and w to z/p, yr, and y^- 

2 



P = 



w 



1 + yp 
2 

l + yr 
2 

1 + Uw 



1, 

1, 

- 1, 



(252) 



and then the final expression for the contribution to the matrix elements from the five-dimensional integral 
is 



{q',l',t',j\M'iA)\q,l,t,j),_^ 



dj3rw 



1 



1 



{l + yrY {l + y^f 
The integral converges nicely in this form 



Icon + IeX + -^IN+EX + luf 



(253) 
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C.2 Miscellaneous Issues 



C.2.1 Rotational Symmetry: j — y —j 

To compute the glucball spcctnim, wc compute the matrix {q' , I' , t' , j\ A4'^{A) \q, I, t, j) and diagonalize 
it for eacli value of j separately. For a given the matrices with j = \j\ and j ~ —\j\ are simply related, 
and we can use this fact to avoid computing and diagonalizing both of them. By inspection, we determine 
that 

{q',l',t',-j\M^{A)\q,l,t,-j}= E il"',l',t',j\M'{A)\q",l,t,j)fiq,q")f{q',q'"), (254) 



where 



fil, q') = ^q,l^q',2 + ^9,2^,1 + ^9,3^,3 " ^g,4^9',4- 



(255) 



This simply means that the basis \q,l,t,—j) is related to the basis \q,l,t,j) by swapping the states 
\l,l,t,j) and \2,l,t,j), and changing the sign of \4:,l,t,j). Renaming basis states and changing their 
phases has no effect on the eigenvalues of the matrix; so |^'(~-'^"(P)) has the same mass as [^'•'"(P)). It 
also means that since 

|*i"(P)) = Y^Ril\q,l,t,j) 

qlt 

= [^m|l'«.*>i)+^'2rj2,/,t,i) + i?^3rtl3,«,t,i) + i?S|4,^,i,i)] , (256) 



we have 



^(-.>(P)^ = E [<t \W,j) + R^, \i,l,tj) + Ril |3,M,j> - RZ |4,M,j> 



J2Ktfil,q')\q',l,t,j). 

q' qlt 



(257) 



Thus by diagonahzing {q' ,1' ,t' ,j\M.'^{A)\q,l,t,j), we also obtain the eigenvalues and eigenstates of 
{q',l',t',-j\M\K)\q,l,t,~3)- 

C.2. 2 Reducing the Number of Matrix Elements to Compute 

There are a few facts that allow us to reduce the number of matrix elements that we have to compute. 
First, because of gluon-exchange symmetry, the basis state \q^ I, t, j) is zero if I + j is even and (/ = 4, or 
ii I + j is odd and 9 7^ 4. Second, A^^(A) is Hermitian, and its matrix elements in this basis are real; so 
it is symmetric in this basis. Third, by inspection, we see that 

{l,l',t',j\M\A)\2,l,t,j)=0. (258) 
Finally, there are some redundancies and additional zeros in the matrix when j = 0: 
{2,l',t',0\M^{A)\2,l,t,0) = {l,l',t',0\M^{A)\l,l,t,0) , 
{l,l',t',0\M^{A)\3,l,t,0) = {2,l',t',0\M^{A)\3,l,t,0) , 
(l,Z',i',0|A4^(A)|4,/,t,0) = -{2,l',t',0\M'^{A)\4,l,t,0), 
(3,Z',i',0|A4^(A)|4,/,t,0) = 0. (259) 
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C.2.2 Estimating Uncertainties in the Spectrum 

When we use the VEGAS Monte Carlo routine to compute the matrix elements, the results for the 
matrix elements have statistical uncertainties. In order to control the resulting uncertainties in the 
spectrum, we would like to have a method that allows us to estimate how accurately we must calculate 
any given matrix element in order for the spectrum to have a desired uncertainty. 

Suppose a diagonal matrix element (g, Z, t,j\M'^{A) \q, I, t,j) is given by 

{qJ,t,j\M\A)\qJ,t,j) ^ Z±S, (260) 

where Z is the Monte Carlo estimate of the matrix element and S is the associated absolute uncertainty. 
Using Mathematica with test matrices, it is straightforward to convince oneself that if 6 is small, then 
it will yield a relative uncertainty e^^^ '-^ S/Z in the eigenvalues of the matrixf^ This translates to a 
relative uncertainty of e^^^ ^ 6/{2Z) in the masses. 

Estimating the uncertainty in the spectrum due to uncertainties in off-diagonal matrix elements is 
more difficult. Using Mathematica with test matrices, the simplest method that we have found that is 
reasonably reliable is to use a type of degenerate perturbation theory. When we have an off-diagonal 
matrix element given by 

{q',l',t',j\M'{A)\q,l,t,j) = Z±5, (261) 

we diagonalize the two matrices 

{q,l,t,j\MHA)\q,l,t,j) Z + 5 

Z + 6 {q',l',t',j\MHA)\q',V,t',j) 

and 

{q,l,t,j\MHA)\q,l,t,j) Z-S 

Z~5 {q',l\t\j\M^iA)\q\l',t',j) 

and we compare their eigenvalues to the eigenvalues of the matrix 
/ {qJ,t,j\M\A)\q,l,tJ) Z 
\ Z {q',l',t',j\MHA)\q',l',t',j) 

We then define e^^^ to be the largest relative deviation that we have found in the eigenvalues, and we 
estimate the resulting relative uncertainty in the mass spectrum to be e^^j^^ ^ e^^^ /2. This estimate tends 
to work well unless there are too many diagonal matrix elements that are nearly degenerate with either 
{q, I, t,j\MHA) |g, I, t, j) or (g', V, t' , j\ M^A) \q' , I', t',j). 

To achieve a relative uncertainty of 0{e) in the glueball masses, we require e^^^ < e for each matrix 
element. This method tends to work reasonably well. The main difficulty is that we have neglected to 
consider the highly nonlinear couplings between the uncertainties in different matrix elements. For this 
reason, as we increase the size of the matrix, the error in our estimate of the uncertainty eventually 
becomes critical. At this point, the spectrum that we get when we diagonalize the matrix becomes 
completely unreliable. (The evidence of the breakdown is sudden contamination of the low-lying wave 
functions with high-order components.) It should be possible to develop more sophisticated methods of 
estimating uncertainties to suppress this problem. 

'^^In the development of this method, we are guided by the principles of quantum-mechanical perturbation theory, although 
we cannot legitimately use perturbation theory to analyze the uncertainties. 
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